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ANALYSIS OF THE DIFFUSE DOMAIN APPROACH FOR A 
BULK-SURFACE COUPLED PDE SYSTEM * 

HELMUT ABELS +, KEI FONG LAM *, AND BJORN STINNER § 

Abstract. We analyse a diffuse interface type approximation, known as the diffuse domain 
approach, of a linear coupled bulk-surface elliptic partial differential system. The well-posedness of 
the diffuse domain approximation is shown using weighted Sobolev spaces and we prove that the 
solution to the diffuse domain approximation converges weakly to the solution of the coupled bulk- 
surface elliptic system as the approximation parameter tends to zero. Moreover, we can show strong 
convergence for the bulk quantity, while for the surface quantity, we can show norm convergence 
and strong convergence in a weighted Sobolev space. Our analysis also covers a second order surface 
elliptic partial differential equation and a bulk elliptic partial differential equation with Dirichlet, 
Neumann and Robin boundary condition. 

Key words, diffuse domain method, weighted Sobolev spaces, well-posedness, diffuse interface 
approximation, asymptotic analysis, bulk-surface elliptic equations. 
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1. Introduction. The diffuse domain approach [39, 31] is a method originating 
from the phase field methodology which approximates partial differential equations 
posed on domains with arbitrary geometries. The method embeds the original do¬ 
main fP c R ra with complicated geometries into a larger domain Q with a simpler 
geometry. Drawing on aspects of the phase field methodology, the diffuse domain 
method replaces the boundary T of fP with an interfacial layer of thickness 0 < e « 1, 
denoted by r e . The original PDEs posed on fP will have to be extended to D and any 
surface quantities or boundary terms on T have to be extended to fields defined on r e . 
The resulting PDE system, which we denote as the diffuse domain approximation, is 
defined on D and will have the same order as the original system defined in IP, but 
with additional terms that approximate the original boundary conditions on T. 

Our model problem will be the following elliptic coupled bulk-surface system: 

-V ■ (AVw) + au = f in IP, 

(CSI) -Vr • (BVrv) + bv + AVu ■ v - [3g on T, 

AVu ■ v = K(v - u) on T. 

Here, VrT and Vr -v denote the surface gradient of v and the surface divergence 
of v on P, respectively. For a precise definition, we refer the reader to [15, Section 
2] or Section 3.2 below. Meanwhile, K, f3 > 0 are non-negative constants, and A = 
( a ij)i<i,j<m and B = ( bij)i<i,j<n denote the matrices with function coefficients a^, 
and bij, respectively. The precise assumptions on the data and the domain will be 
given in Section 2.1. 

We now embed O* uT into a larger domain D c R n . The location of the original 
boundary P is encoded in an order parameter tp as its zero-level set. A typical choice 
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for ip is a function of the signed distance function of T, d : Cl -»■ M, which is defined as 


( 1 . 1 ) 


d(x) 


-inf zer |a; - z\ 
■ 0 

inf zer \x - z\ 


for x e Cl*, 
for itf, 
for x e Cl \ Cl*. 


Let x'fi* denote the characteristic function of Cl* and Sr '■= 'H n_1 l T denote the 
Hausdorff measure restricted to T (see Theorem A.l for a definition in the sense of 
distributions given by the Dirac measures of Cl* and T). An equivalent distribution 
form for (CSI) is (see Appendix, Theorem A.l) 

-V • (xn*AVu) + xn*au = xn*f + S r K(v - u) in 1 ?'(fl), 

-V • ( SrBVv) + Srbv = Srfdg - SrK(y - u ) in V' (Cl). 

The diffuse domain approximation of (CSI) is derived by approximating xn * and 
Sr with more regular functions £ E (<^), 5 s (tp), indexed by e, which is related to the 
thickness of the interfacial layer T e . In order words, the diffuse domain approximation 
of (CSI) is 

-V • ( £ E A E '\7u e ) + £, e a E u e - £ e f E + S e K(v e - u e ) in 0, 

(CDD) ^ ( 5 e B E X/v e ) + 5 e b E v e = S e /3g E - 5 e K(v e - u e ) in D, 


where A E ,B E ,a E ,b E , f E and g E denote suitable extensions of A,B,a,b,f and g to 
the larger domain Cl. The precise assumptions on the extensions will be outlined in 
Section 2.4. 

Formally, £ e (<£>) -»■ xo*, S E (p) -»■ <5r as e -* 0, and so in the limit of vanishing 
interfacial thickness, we recover the distributional form for (CSI). 

For here onwards, we refer to the original problem posed on Cl* as the SI (sharp 
interface) problem, and the corresponding diffuse domain approximation posed on Cl 
will be denoted as the DD (diffuse domain) problem. 

We remark that replacing the original boundary T with an interfacial layer T e 
transforms the diffuse domain approximation into a two-scale problem. This is sim¬ 
ilar to phase field approximations of free boundary problems. The idea of adopting 
the phase field methodology to approximate partial differential equations has been ap¬ 
plied to study diffusion inside a stationary cell [26], Turning patterns on a membrane 
[30], wave propagation in the heart [12, 20], two-phase flow [3] and soluble surfac¬ 
tants [22, 42]. In [9, 10], a diffuse domain type method, denoted as the smoothed 
boundary method, has been applied to solve partial differential equations on irregular 
domains with homogeneous Neumann boundary conditions using spectral methods. 
A generalised formulation of the smoothed boundary method can be found in [44]. 

The phase field methodology provides us with two candidates for p(x). The first 
is based on the smooth double-well potential ipDwiv) = ^(l - ^ 2 ) 2 and leads to 


Pdw(x) : = tanh 



for x eCl. 


The other is based on the double-obstacle potential [6, 13] 


i’Doiv) = ^ “ ^ 2 ) + h- i,i](^)> 7 [-i,i](^) 



if [-1,1], 
otherwise, 




ANALYSIS OF THE DIFFUSE DOMAIN APPROACH 


3 


and leads to 


Vdo(x) 


+ 1 , 

- sin(d(a;)/£), 

- 1 , 


if d(x ) > ef, 
if \d(x)\ < e§, 
if d(x) < -ef • 


Note that in the case of the double-obstacle potential, the interfacial layer T e has finite 
thickness of £7r. We remark that for general phase field models ipDW or ip do is the 
leading order approximation for the order parameter ip. In our setting, the location of 
the boundary T is known and hence we can use ipDW or ipDO as the order parameter. 
A common regularisation of xn* based on the smooth double-well potential used in 
[31, 41, 42] is 


d 1} (*) = \o- - <pdw(x)) = ^ |i - tanh (f/f;))' 

While an alternative based on the double-obstacle potential is 

£ e 2 ) 0 ) = ^(1 - Pdo(x)). 


There are many regularisations of (5r available from the literature [42, 17, 39, 28]. One 
well known approximation of 5r is a multiple of the Ginzburg- Landau energy density 
(see [37]) 


| Iv^l 2 + 

2 e 

From the above discussions regarding the double-well and the double-obstacle poten¬ 
tials, we have two candidates for the regularisation to <5r: 


4 ' ,(l) -i7S 8ed,4 (75)' eHx) 


^ 2 
-COS 

7re 



X{xzQ : |d(x)|<e-|-} • 


The convergence analysis of the diffuse domain approach (with the smooth double¬ 
well potential), in the limit e -* 0, has only been done in the context of recovering 
the original equations via formally matched asymptotics (see [31, 41, 42, 29]). A 
first analytical treatment of convergence in one dimension and on a half-plane in two 
dimension can be found in [21], where the error between the solution to a second order 
system and the diffuse domain approximation in the L°° norm is of order 0(e 1_At ), 
where p > 0 arbitrarily small. 

In [16], a diffuse domain type approximation for an advection-diffusion equation 
posed on evolving surfaces is considered. Motivated by modelling and numerical sim¬ 
ulations, the diffuse domain approximation utilises a double-obstacle type regularisa¬ 
tion. Note that the regularisations from the double-obstacle potential are degenerate 
in certain parts of the larger domain f2 (in particular they are zero outside T e for 
the set-up in [16]). Consequently the corresponding diffuse domain approximation 
becomes a degenerate equation and weighted Sobolev spaces are employed. The chief 
results in [16] are the well-posedness of the diffuse domain approximation and weak 
convergence to the solution of the original system. 

For a bulk second order elliptic boundary value problem, the convergence analysis 
has been studied in [11]. There, a double-obstacle type diffuse domain approximation 
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to a Robin boundary value problem is studied, with the specific choice S £ = |V£ e |. 
The authors are able to deduce trace theorems, embedding theorems and Poincare 
inequalities for Sobolev spaces weighted with Under suitable assumptions on the 
decay of £ e near <9T £ , the authors showed that there exists a p > 2 such that the error 
between the solution to the diffuse domain approximation and the solution to the 
Robin problem is of order 0(e 5_ p) in a weighted Sobolev norm. Similar results for 
the Neumann and Dirichlet boundary conditions are also given. 

In this work, we show that weak solutions to (CDD) converge to the unique weak 
solution to (CSI) under appropriate assumptions. Thanks to the property of the 
problem, we are also able to show strong convergence in H 1 ( Cl*) for the bulk quantity, 
while for the surface quantity we have norm convergence and strong convergence in 
a weighted Sobolev space. The techniques we develop here can also be applied to a 
second order surface linear elliptic equation and a second order linear elliptic equation 
with Robin, Neumann and Dirichlet boundary conditions. 

The structure of this article is as follows: In Section 2 we introduce the assump¬ 
tions and the main results. In Section 3 we prove several technical results that will 
simplify the proof of the main results, which are contained in Section 4. In Section 5, 
we compare our results with the results in [16, 11] and discuss some insights and issues 
we encountered in the analysis. In the appendix we collect a few results regarding the 
distributional formulations and properties of the trace operator. 

2. General assumptions and main results. We observe that setting A = 
0 nx ™, a = 0, / = 0, K = 0, /3 = 1 in (CSI) leads to an elliptic equation on T, which we 
denote as (SSI): 


(SSI) - Vr • (BVrv) + bv = g on T. 


Meanwhile, setting B = 0”* ra , b(x) = (3 and formally sending K -*■ oo in (CSI) leads to 
the Robin boundary condition (RSI): 


(RSI) 


-V • (AVu) + au = f in Cl, 
AVu ■ v + f3u = pg on T. 


Neglecting the second equation in (CSI), setting v = g and formally sending K -*■ oo 
gives the Dirichlet boundary condition (DSI): 


(DSI) 


-V ■ (A\7u) + au = f in SI, 
u = g on T, 


while the Neumann boundary condition (NSI) can be obtained by setting B = 0 nxn , 
b = 0, (3 = 1 and neglecting the last boundary condition in (CSI): 


(NSI) 


-V ■ (*4Vu) + au = f in SI, 
AVu -v = g on T. 


In order to cover these derived sharp interface problems we make the following as¬ 
sumptions. 

2.1. Assumptions on the data. 

Assumption 2.1 (Assumptions on domain). We assume that Cl* is an open 
bounded domain in R” with compact C 2 boundary T and outward unit normal v. Let 
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f l be an open bounded domain in R n with Lipschitz boundary dfl such that Q* c Q 
and r n dfl = 0 . In addition, for (CSI) and (SSI), we assume that T is of class C 3 . 
Assumption 2.2 (Assumptions on the data). We assume that for 1 < i, j <n, 

ciij,a e /eL 2 (fT), b tj ,b e L°°(r), geL 2 (r), 

and there exist positive constants 9 q, 61,82, 03 such that 

(A(x)Ci)-Ci>flo|Ci| 2 , (B(p) C2) • C2 > 6>1 IC2I 2 , a(x)>6 2 , b(p) > 0 3 , 

for all x € H*, p € r, Ci e R" and C 2 e T P T c R n . 

Assumption 2.3 (Specific assumption for (DSI)). For (DSI), we assume 

5 eUi(r). 


Here, for a bounded Lipschitz domain D c R”, the fractional-order Sobolev spaces 
H?{dD) is defined as 


H^(dD):=\ueL 2 (dD): f f 

{ JdD JdD 


\u(x)-u(y)\ 


■y\ r ' 


dH"- 1 ( 1 ) dH"' 1 (y) < < 


equipped with the inner product 


(u,v) 


H?(dD) J dD 


■-L 


uv dlC 


[ [ 

JdD JdD 


(u( x)-u(y))(v(x)-v(y)) 
I x-y\ n 


dU n - 1 (x)dU n ~ l (y), 


where d'H ” -1 denotes the n— 1 dimensional Hausdroff measure. 

Assumption 2.4 (Specific assumption for (NSI)). For (NSI), if g * 0 then we 
assume that 


g e H 2 (r), and a ^ e C° :1 (H*) for 1 < i,j < n. 


2.2. Results on the sharp interface (SI) problems. Let 70 : W 1,1 ( S2*) 
L 1 (r) denote the trace operator. If T is Lipschitz and g e Hz(T), then by Theorem 
B.l, there exists g e 7L 1 (fT‘) such that 70 (g) = g , and (DSI) is equivalent to the 
following homogeneous problem (DSIH): 


(DSIH) 


-V • (AVw) + aw = / + V • (AVg) - ag in H *, 
w = 0 on r. 


If w is a weak solution to (DSIH), then u := w + g is a weak solution to (DSI). 

We can also consider (DSIH) as the limiting problem /3 -* 00 of the following 
homogeneous Robin problem (RSIH): 


(RSIH) 


-V • (AVw) + aw = / + V • (AVg) - ag in f2*, 
AVw ■ v + fiw = 0 on T. 


In fact, we have (see also [32]): 
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Lemma 2.1. Let Q* be an open bounded domain with C 1 boundary T. Let the 
data satisfy Assumption 2.2 with g e For each /3 > 0, (RSIH) is well-posed 

with unique weak solution w^ € H 1 (fl*), and (DSIH) is well-posed with unique weak 
solution wd e (see Theorem 2.1 below for the well-posedness). Then, as 

/3 -»■ oo, converges weakly to w B in H 1 ^ 2*). 

The proof of Lemma 2.1 will be outlined in Section 4.1. 

Similarly, by Assumption 2.4 and Theorem B.2, there exists h e H 2 (il*) such 
that AVh■ v - g in the sense of traces. Then, (NSI) is equivalent to the following 
homogeneous problem (NSIH): 


(NSIH) 


-V • (AVw) + aw = / + V • (AVh) - ah in f 1*, 
AVw ■ v - 0 on T. 


If w is a weak solution to (NSIH), then u := w + h is a weak solution to (NSI). We 
point out that, if h t 0, then the coefficient A in (NSIH) is as assumed in Assumption 
2.4, i.e., aij e C 0,1 (f2*) for 1 < i,j < n. 

We now state the well-posedness results for the SI problems. For convenience, 
let us define the following bilinear forms as '■ H 1 (I1*) x H 1 (Q*) -*■ R, as ■ iL^T) x 
7L 1 (r) K, l B : L 2 (fT) x L 2 (fl*) R, and l s ■ L 2 (T) x L 2 (T) R: 

(2.1) a B ((p,ip) := / AX/ip-yip + atpip dx, l B (ip,ip) ■= / ifiip dx, 

J n* Jn* 

(2.2) a s (<p,ip) := Vrip + b(pipdTL n ~ 1 , l s (ip,ip) := J^ipipdH^ 1 . 


Theorem 2.1 (Well-posedness for the SI problems). Let Assumptions 2.1, 2.2, 
2.3 and 2-4 be satisfied. Then there exist unique weak solutions 


( u, v ) e H\ fT) x H\T) for (CSI), 
u R eH\n*) for (RSI), 
u N €H\n*) for (NSI), 


v s eH\T) for (SSI), 
w R eH\n*) for (RSIH), 
WncH 1 ^*) for (NSIH), 
w D eH((Ll*) for (DSIH), 


such that for all {p e 7L 1 (S 1*), ipo e ip e 7L 1 (T), 


a B {u,tp) + a s (v, ip) + Kl s (v - 70 (u),ip- 7 o(<^)) = h(f,T) + Ph(g,tp), 

as(v s ,tp ) = ls(g,ip ), 

a B (u R ,ip ) +^/s( 7 o(u fl ), 7 o(<^)) = l B (f,<p) + /%s(5,7o(^)), 
a B (w R , p) + /?/s(7o(uTi),7o(^)) = h(f, p) ~ a B (g, ip), 

a B (u N ,p) =l B (f,<p) + l S (g, Ao(t)), 
a B (w N ,p ) = l B (f+y -(Ayh)-ah,ip), 
a B (w D ,ip 0 ) = l B (f,p 0 ) ~a B (g,p 0 ). 


Moreover, there exist constants C, independent of(u,v), vs, ur, wr, w b , un, wn, 
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such that 


II u IIh 1 (q*) + IMIn^r) ^ C(ll/ili2(n*) + llffli 2 (r)); 

ll w s , ||fl'i(r) ^ Cllflllz, 2 (rp 

II^RlfPCn*) + ^IIto('^-r) llz, 2 ( r ) - ^'(\\f\\ 2 L 2 (u*) + II^IIz, 2 (r))^ 

Il u ’-Rlli7 1 (fi*) + Pi'yoi'WR.) IIL 2 (r) - C(ll/ll l 2 (n*) + II3II h 1 (n*))’ 
ll u Y||^i(n») - C(l/lll 2 (n») + llfl , lll 2 (r))> 
ll' u dvlliri(n*) ^ C(l/lll 2 (n») + ll^llH 2 (n*)); 

ll w D||^i(n*) - C(ll/lli 2 (n») + II3II Ip (n*))- 

We note that the constant C for wn depends on ||.A|| c .o,i(fj*p while the constant 
for w R is independent of f3 (see Section 4.1). 

We will outline the proof for (CSI) in Section 4.2. The well-posedness for the 
other SI problems follow similarly from the application of the Lax-Milgram theorem, 
and hence we will omit the details. 

2.3. Derivation of the diffuse domain approximations. Motivated by The¬ 
orem A.l, replacing xci* and Sr by £ e and S e , we obtain the corresponding diffuse 
domain approximations for (SSI), (RSI), (RSIH), (NSI) and (NSIH), which we denote 


by (SDD), (RDD), (RDDH), (NDD) and (NDDH): 

(SDD) - V • (5 e B E Vv e s ) + S e b E v e s = S e g E in SI, 

(RDD) - V • (£ e A E \7u R ) + £ E a E u R + /35 e u R = £ e f E + fdS e g E in 0, 

(RDDH) - V • (£ e M £ V«4) + £ e a E w £ R + pS e w e R = $ e {f E - a E g E ) + V • ^ e A E \/~g E ) in SI, 

(NDD) - V • (£ e M B Vi4) + Cea E u e N = f e f E + S e g E in SI, 

(NDDH) - V • (£ e A e \7w e n ) + f, e a E w e N = f, e f E + &V ■ ( A E Vh E ) - f, e a E h E in SI. 

Here, terms with the superscript E denote extensions to the larger domain D. We 
will outline the specific choices for these extensions in Section 2.4 below. 


We note that Theorem A.l does not allow us to infer an approximation for PDEs 
with Dirichlet boundary conditions. However, motivated by Lemma 2.1, we obtain 
an approximation by scaling /3 appropriately in (RDDH). Hence, a diffuse domain 
approximation to (DSIH) is (DDDH): 

• (ZeA E X7w e d ) + £ e a E w e D + -S E w e D = &(/ B - a E g E ) + V • (&A E Vg E ) in SI. 

£ 

Here, one may infer that if w e D converges to some function w° as e -*■ 0, then S e w e D 
would converge to the trace of w°. So the term ^S e w e D can be seen as a penalisation 
to enforce 70 (re 0 ) =0 in the limit e -*■ 0 . 

Remark 2.1. We could set /3 = — f or m > 0 to obtain a family of diffuse domain 
approximations to (DSIH). Then (DDDH) is the case m = 1. We refer the reader to 
the discussion in Section 5.1 for more details. 

Remark 2.2. It will turn out that the problem (NDDH) has better properties than 
(NDD). In particular, the subsequent analysis for the Neumann problem is performed 
with (NDDH) instead o/(NDD). We refer the reader to the discussion in Section 5.2 
for more details. 
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Finally, we assume a homogeneous Neumann boundary condition for our diffuse 
domain approximations: 

A E X/u e ■ VQQ. = A E Vw e ■ Von = B e X/v e ■ vqq. = 0 on dfl. 

2.4. Assumptions on data extensions. In general, the extension operator is 
not unique. For Sobolev functions / e W l,p {i 1*), where l > 1,1 < p < oo, one can use 
an order l - 1 reflection about the boundary T to extend / to the exterior of f l*. The 
case l = 1 is outlined in the extension theorem [18, Theorem 1, p. 254], while for 
higher order reflections, we refer to [2, Theorem 5.19, p. 148] or [34, p. 43-44]. 

Thus, for (DSIH) we can extend g e H 1 ^*) to g E e 14 1 (fl) by the extension 
theorem, and for (NSIH) we can extend h e ff 2 (H*) to h E e H 2 (tt) by the reflection 
method. 

For the other bulk data, i.e., A, a and /, we make the following assumption: 

Assumption 2.5 (Extension of bulk data). Let 1 < i,j < n. For aij 7 a e L°°(fl*) 
and f e L 2 (Q*), we assume that there exist extensions a E 7 a E e L°°(Q), and f E € 
L 2 (Ll) such that A e = ( a E ) l <i,j<n is uniformly elliptic with constant 9 q and a E (x ) > 62 
for a.e. x e Cl. 

For (NSIH) we additionally assume that the extension a E € C'°’ 1 (fl) for 1 < i, j < n 
and the matrix A E preserves uniform ellipticity. 

The surface data will be extended in a specific way. Recall that by Assumption 
2.1, r is at least of class C 2 . We define the tubular neighbourhood Tub r (T) of T with 
width r > 0 as 


Tub r (r) := (x e f!: |d(a;)| < r). 

Then, by [24, Lemma 14.16], there exists 77 > 0 such that the signed distance function 
d to r is of class C , 2 (Tub I? (r)). Moreover, it can be shown that d is globally Lipschitz 
with constant 1 (see [24, Section 14.6]). 

For each y e T, let T y T and v(y) denote its tangent space and outward pointing 
unit normal, respectively. A standard result in differential geometry shows that for 
77 sufficiently small, there is a C 2 diffeomorphism between Tub r, (T) and T x (- 77 , 17 ) 
given by 

(2.3) 0 ":Tub"(r)-rx (-77,77), Q v (x) = (p(x),d(x)), 

where, for any x e Tub ?? (r), we define the closest point operator (see [36] or [15, 
Lemma 2.8]) p : Tub 7 ' 1 (F) T by 

(2.4) p{x) := x - d(x)v{p(x)). 

Then, we also have 

(2.5) Vd(x) = v(p(x)) for x € Tub r, (r). 

Definition 2.1 (Constant extension in the normal direction). For any g e L q (T), 
1 < q < 00 , we define its constant extension g e ojfT to Tub ?, (r) in the normal direction 
as 

(2.6) g e ( x ) - g(p( x )) f or x € Tub r, (r). 

By Corollary 3.1 below, we have that g e e L 9 (Tub r, (r)) if g e L^T) for 1 < q < 00 . 
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For the surface data, we make the following assumption: 

Assumption 2.6 (Extension of surface data). Let 1 < i.j < n. For bij,b e L°°{ T) 
and g £ L 2 {T ), let bT,b e £ L°°(Tub T, (r)) and g e £ L 2 (Tub ?J (r)) denote the constant 
extensions ofbij,b and g offT to Tub ,; (r) in the normal direction, respectively. We 
assume that there exist extensions b E ,b E £ L°°(Tl), and g E £ L 2 (Tl) of bfj,b e and 
g e , respectively, such that B E = (b E )i<ij< n is uniformly elliptic with constant B\ and 
b E (x ) > 03 for a.e. x £ ft. 

Remark 2.3. To avoid confusion with the different types of extensions for the 
bulk and surface data introduced above, we will make the following clarification and 
use the notation: The superscripts Ea, Er and Ec will denote arbitrary extension, 
extension by reflection and constant extension in the normal direction. More precisely, 

• For 1 < i,j < n, the extensions b Ec ,b Ec £ that appear in (CDD) 

and (SDD) are the extensions ofbij,b £ L°°(r) constructed from extending 
constantly in the normal direction offT to Tub r, (r) ; and then arbitrary from 
Tub ?, (r) to fl, such that Assumption 2.6 is satisfied. 

• The extensions g Er £ H 1 ^) and h Er e H 2 (Tl) that appear in (DDDH) and 
(NDDH), respectively, are extensions by reflection of g £ Ff 1 ^*) and h £ 
H 2 {TL*). 

• The extension g Ec £ L 2 (TL) that appears in (CDD), (SDD) and (RDD) is the 
extension of g £ L 2 (T) constructed from extending constantly in the normal 
direction offT to Tub’ , (T), and then arbitrary from Tub ?, (r) to Q. 

• The extensions a Ea £ and f Ea £ E 2 (Tl) that appear in (CDD), (RDD), 

(NDDH) and (DDDH) are arbitrary extensions of a £ and f £ 

L 2 (fl*) such that a Ea > O 2 a.e. in fl. 

• For 1 < i,j < n, the extension a Ea £ L°°(fl) that appears in (CDD), (RDD) 

and (DDDH) is an arbitrary extension of aij £ such that the matrix 

A Ea preserves uniform ellipticity. Meanwhile, the extension a Er £ C 0 , 1 (H) 
that appears in (NDDH) is the extension by reflection of £ C' 0 , 1 (H*), and 
the matrix A Er preserves uniform ellipticity. 

2.5. Assumptions on regularisations. We first introduce the functions £ and 
S, from which the regularisations and S e are constructed by a rescaling. 

Assumption 2.7. We assume that f : R -»■ [0,1] is a continuous, nonnegative, 
monotone function such that 

if x < 0 , 
if x > 0, 
if x = 0 . 


(2.7) 0 < £(t) < £(0) = - < £(s) < 1, Vs<0 <t, and lim^( —) 

2 e-»o \ e ) 


Assumption 2.8. We assume that 6 '■ R -»■ R>o is a C 1 , nonnegative, even 
function such that 

(2.8) [ S(s) ds = 1 , 6 (si) > 6(s 2 ) if |si| < |s 2 |, 

Jr 

(2.9) [ ^Lds+ f y/S{a)+S(a)(\s\ + \s\ 2 )ds =: C 5 , int < 00 , 

J{sgR:<5(s)>0} (5(s) Jr 

and for any q > 1 , 


lim —S 

e—* o e q 



( 2 . 10 ) 


if x + 0 , 
if x = 0 . 
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Moreover, we assume there exists a constant > 0 such that 
(2.11) C^S(t) < f(t) for all t e R. 


Definition 2.2. Let d(x) denote the signed distance function to T. For x € Cl 
and for each e € (0,1], we define 


( 2 . 12 ) 

with 

(2.13) 


£e(z) : =£ 




Cl e := {x e Cl : f E (x) > 0}, r e := {x e it: S e (x) > 0}. 


By (2.7), (2.10) and (2.11), we observe that 

Cl* u r c fl £ , r c T £ c Cl e , for all £ > 0. 


One can check that Assumptions 2.7 and 2.8 are satisfied by our candidate regu- 
larisations originating from the double-well potential: 



and from the double-obstacle potential: 

i D °{x) := X(-oo,-§ )(x) + i(l - sin(x))x[-f ,f](aO, 

(2-14) 2 

S DO (x) := -cos 2 (a;)x[- f , f ](a:). 

Remark 2.4. For any e > 0, and 6^ w derived from the double-well potential 
are non-degenerate in Cl, i.e., Cl e = T e = Cl for all e > 0. However, and 5®° 
originating from the double-obstacle potential are degenerate in Cl. In particular, for 
(2.14), 

a e = fi*uTub £ 5 (r), r £ = Tub £ - 5 (r). 


Moreover, if £\ < £ 2 , we have fl £l c Cl e2 and T £l c T £2 . However, the framework 
of weighted Sobolev spaces is flexible enough to allow us to deduce well-posedness 
of the diffuse domain approximations with both the double-well and double-obstacle 
regularisations. 

2.6. Weighted Sobolev spaces. Due to the presence of and 5 e in the diffuse 
domain approximations, the natural function spaces to look for well-posedness are 
Sobolev spaces weighted by and S e . In the following, measurability and almost 
everywhere are with respect to the Lebesgue measure. 

Definition 2.3. For fixed £> 0, we define 


L 2 (ti E ,t; s )-.= 

i 2 (r e ,4):= 


f:Cl e - 

R measurable s.t. 

fa, (, l/l 2 dx < 00J 

/ = r e - 

* R measurable s.t. 

J r Se I /] 2 dx < 00 j 
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By Lipschitz continuity of the signed distance function d, and the continuity of 
£(•), we see that £ E (-) is continuous, and consequently j- is bounded in all compact 
sets B c S\. Thus, j- € Lj oc (Q e ) and by Holder’s inequality we have the continuous 
embedding L 2 (fl E ,t; E ) c Lf oc (tt e ) (see also [27, Theorem 1.5]). Thus, we can define 
derivatives for / e L 2 (fl e ,£ £ ) in a distributional sense. I.e., for any multiindex a, we 
call a function g the a th distributional derivative of /, and write g = D a f, if for every 


f fD a <pdx = (-l)H f get *d 
j£i e Jq £ 


x. 


We define the vector space 


W 1 ’ 2 ^,^) := {/ e L 2 (D E ,£ E ) : D a f e L 2 (D E ,£ E ) for |a| = 1}. 


A similar definition for the vector space W 1 ’ 2 (T e ,<5 e ) can be made since J- e 
L L c ( r £) b y the continuity of 5 e {-). 

For the subsequent analysis, we will present the proofs with the double-well reg- 
ularisation in mind and detail any necessary modifications for the double-obstacle 
regularisation afterwards. 

To streamline the presentation, it is more convenient to have a fixed domain when 
working with weighted Sobolev spaces, hence we introduce the following notation: 

Definition 2.4. For fixed e > 0, we define 

L 2 { S2,£ e ) := {/ : D -»■ R measurable s.t. f |n E e T 2 (fl E ,£ E )}, 

H 1 { f2,£ e ) := {/ : f2 -»■ R. measurable s.t. f |o e c W 1 ,2 (fi e ,£ e )}, 
with inner product and induced norms: 


(f,g)L 2 ( n,£ e ) : - 

f f.efg dx = f £ e fg dx, 

J n Jn e 

/!o,|e (fif)L 2 (n,£ e ): 

(/)S , )i/ 1 (n,£ e ) : = 

[ Ze(fg + Vf-Vg)dx, 

JSl 

ll/I i,j e := (/>/)j? 1 (n,£ e ) 

: use the identification: 



f = g <=> /(x) = 5 (x) a.e. 

in fl E . 


A similar notation and identification are used for L 2 (fl,S E ) and H 1 (H,6 e ). Further¬ 
more, we define 

V e := {/ : fl -»■ R measurable s.t. f |n e e hF 1 , 2 (fl £ , £ e )} with 
(f,g)v e : = f n ts(fg + Vf-Vg) + Sefgdx, 

and 


W E := {/ : fi -*■ M measurable s.t. f |n e e TT 1,2 (H e , £ e )} with 

(f,g)we : = f Ce(/3 + v/- Vg) + -S E fgdx. 

Jn e 


Remark 2.5. We observe that i/ 1 (f2,£ E ), V E and W E are the same vector spaces 
but with different inner products and induced norms. It will turn out that the well- 
posedness to (RDD), (DDDH) and (NDDH) depend critically on the choice of the 
inner product. 
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2.7. Results on the diffuse domain (DD) problems. Similar to the above, 
we introduce the following bilinear forms a e B : x R, a% : 

iJ 1 (fl,5 £ ) x 17 1 (f l,5 e ) -*■ R, l e B : L 2 (H,£ £ ) x L 2 (H,£ £ ) -»■ K, and l E s : L 2 (Ll,8 e ) x 
L 2 (n,<5 £ ) -»■ R: 

(2.15) a B ((p,ip) := f £, e (A Ea \7(f-Vip + a Ea <pip) dx, l e B (tp,ip) := [ £ e (pip dx, 

v/ Q w/ Q 

(2.16) a E s (<p,ip) := f 6 e (B Ec \7<p ■ X/ip + b Ec <pip) dx, l e s (ip,ip) := f 8 e ipip dx. 

v/ $T2 O 

For (NDDH), we replace in afj( v ) with „4 Br . 

Our first main result is the well-posedness of the diffuse domain approximations: 
Theorem 2.2 (Well-posedness for the DD problems). Let tt c R n be an open 
bounded domain with Lipschitz boundary. Suppose Assumptions 2.5, 2.6, 2.7 and 2.8 
are satisfied. In addition, for (CDD) we assume that 0$ > K. Then, for each e > 0, 
there exist unique weak solutions 

(u E ,v E ) eV E xH\n, S E ) for (CDD), «| e H 1 ^^) for (SDD) 

u e R e V £ for (RDD), w E N e H\ fi,&) for (NDDH), 
w e D e W £ for (DDDH), 

such that for all ip e R 1 (D, £ £ ), ip e 17 1 (D, 5 e ), (p e V £ , $ e W £ , 

(2.17) a E B (u E , <p) + a E s (v e ,ip) + Kl E s (v E -u E ,ip-<p) = l E B (f Ea ,<p ) + /3l E s (g Bc ,ip), 

(2.18) a E s (v E s , ip) = l e g{g Ec , ip), 

(2.19) a E B (u E R , <p) + f3l E s (u R , <p ) = l B (f Ea , (p) + fil e s (g Ec , <p), 

( 2 . 20 ) <&(«& + $) + ^|(^, $) = $), 

(2.21) a E B (w E N , ip) - Z|(V • (^V/^) - p) = l E B (f Ea , p). 


Moreover, the weak solutions satisfy 


( 2 . 22 ) 


1.9 L 2 (r)), 

(2.23) 

U U !I \l,S e - C \Ifli L 2 (r)> 


(2.24) 


9 L 2 (r))> 

(2.25) 

\\ w d\\i,^ + Wd\ 2 q 1 ^ C{\\f Ea \\ 2 L 2 {n) + 

\\~9 Er \W) 

(2.26) 

\\™%\\h £ <c(\\f Ea \\ 2 LHn) + 

llOl^n) 


where the constants C are independent of e. 

We point out that, thanks to the fact that g Ec is the constant extension of g in 
the normal direction, the estimates in (2.22), (2.23) and (2.24) are independent of e. 

2 .8. Convergence results. Recall that the unique weak solutions to (CS1), 
(SSI), (RSI), (DSIH) and (NSIH) are denoted by (u,v), vs, u R , wd and wn, re¬ 
spectively. Keeping the notation that the unique weak solutions to (CDD), (SDD), 
(RDD), (DDDH) and (NDDH) are denoted by (u e ,v e ), v E s , u R , w e D , and w E N , our 
second main result in the weak convergence of diffuse domain approximations: 

Theorem 2.3 (Weak convergence of diffuse domain approximation). Suppose 
Assumptions 2.1, 2.2, 2.3, 2.4, 2.5, 2.6, 2.7 , and 2.8 are satisfied. In addition, for 
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(CDD), we assume 0 3 > K. Then, as e 0, 

u e |n», u e R |n*, w e D |n*, w% |n* converge weakly to u, ur, wd, % m H 1 (Cl*), 

respectively. While, for any tf> e iJ" 1 (r) with extension iJ j Ec € i/ 1 (f2) (as constructed 
in Corollary 3.1 below), as e -*■ 0, 

f S e v £ i/j Ec dx^ f vipd'H n ~ 1 , f S E Vv E ■ Vip Ec dx -*• f Vrv ■ Vrf’d'H ^ 1 , 

J n Jr Jn Jr 

J' SeVri/j^ dx-> J vs’4>d’H n ~ 1 , J S E \/v e s ■ Vif Ec dx^ J S/rvs • Vr^d"H" _1 . 
Let 

(2.27) Cr := {ip Ec e iL 1 (f2) : if Ec constructed from if e JJ 1 (r) as in Corollary 3.1} . 

Then, one can formally interpret the above weak convergence of the surface quantities 
in the following way: 

5 e v E -- 5ru, 8 e Vv £ -- drVr'c in Cf as e 0. 

The precise definition of the above notion of weak convergence for the surface quan¬ 
tities can be found in the proof of Lemma 3.9. In fact, in Lemma 3.9, the weak 
convergence for the lower order term holds for all iJ 1 (fl) functions: As £ -*■ 0, 

f 5 e v £ (pdx -* f v r )o(g))d'H n ~ 1 for all tp e H 1 (ft). 

Jn Jr 

Thanks to the coercivity of the bilinear forms a e B and a E s , we also obtain strong 
convergence: 

Theorem 2.4 (Strong convergence of diffuse domain approximations). Suppose 
Assumptions 2.1, 2.2, 2.3, 2.4, 2.5, 2.6, 2.7, and 2.8 are satisfied. In addition, for 
(CDD) we assume 0 3 > K. 

Let v Ec ,v Ec e U 1 ( fl) denote the constant extensions of v,vs g U 1 (T) in the 
normal direction, as constructed in Corollary 3.1. Let u Er ,uf> r e iL 1 (f2) denote the 
extensions of u,ur g H 1 (D >1 ‘) by reflection. Then, as e ^ 0, 

IK In* -u||ifi(n*) + IK ~ uEr lo,s e + IK ~ v Ec \i,6e 0’ 

|v!-«f c lk«. -*- 0 , 

II u r In* -Mfl||ffi(n*) + II Ur - u R r ||o,4 -»■ 0, 

\\w e N In* -WAr|ffi(n*) 0, 

\\w e D |n* -wdWh 1 ^*) -*■ 0. 

Consequently, we also have the norm convergence: As e -*■ 0, 

|K||o,5 e -»■ |7o(«)IK(r), ||« e ||i, 4 . -* IMK(r), 
ll u slli, Se IKIK(r)> Kllo,«. -*■ Ito(hr)|| z, 2 (r) ■ 

In the next section, we present some technical results that will be essential to 
prove the main results. 
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3. Technical results. 

3.1. Change of variables in the tubular neigbourhood. Choose 77 > 0 and 

the diffeomorphism 0 ,f : Tub r, (T) -+Tx (- 77 , 77 ) defined in (2.3). For t e (- 77 , 77 ), let T t 
denote the level set {x e : d(x) = t). Then, by the co-area formula [18, Theorem 5, 
p. 629] or [19, Theorem 2, p. 117], and that ]Vd| = \i/\ = 1, we can write 

(3.1) f f{x) |Vd(x)| dx = f f(x) dx = f f /d'H” _1 dt. 
JTub’i(r) v v ' JTubi(r) v ' J~ v Jv t 

We define the mapping p t : T -*■ T t by 

(3.2) p t (p) = p + tu(p) for peT. 

This map is well-defined and is injective due to the diffeomorphism Q 11 . Then, by a 
change of variables, we obtain 

f fdH "- 1 = f f(p + Mp))\det((Vp t ) T (Vp t ))fi dTT" 1 , 

Jr t Jr 

where V pt is the Jacobian matrix of p t . To identify det((Vpt) T (Vpt)) as a function 
of t we use local coordinates. 

Since T is a compact hypersurface, we can always find a finite open cover of T 
consisting of open sets Wi c R", 1 < i < N such that T c (jf^i Wj. For each 1 < i < N, 
let a.i(s ) denote a regular parameterisation of Wi nr with parameter domain «S ? cM n_1 , 
i.e, a.i : Si -»■ Wi n T is a local regular parameterisation of T (the existence of such 
local regular parameterisations follows from the regularity of T). Let 

Ji,q(s) := (5 Sl ai(s),... ,d Sn _ 1 a i (s),7x(a i (s)) e R n!m . 

Since at is a regular parameterisation, the tangent vectors {d Sj ai}i<j< n -i are linearly 
independent and hence the determinant det Jtp * 0. Then, for any / e L) oc (r), 

(3.3) f fdU n ~ x = f f(a{(s)) |det Ji j0 | ds. 

JWinT J Si 

By the injectivity of pt, Ft is also a compact hypersurface with a finite open cover 
{pt(Wi n r)}^. In addition, p t ° cm is a local parameterisation of r t . Let 

Ji, v (s,t ) := (d Sl at(s) + td Sl i'(a i (s)),.. .,d Sn _ 1 a i (s) + td Sn _ 1 iy(a i (s)), v(a t (s)) e R nxn 
= Ji,o(s) + t(d Sl v(on(s )),... ,a s „_ 1 jx(a i (s)),0) =: A,o(s) + tBi(s). 


A short calculation shows that 

det J itV = det( Ji j0 + tBi) = (det J it 0 )(det(I + = (det Jip)t n det (f -1 / + J^pBi) 

= (det J ii0 )t n (i - -^-tr (-J-IB,) + ■■■ + (-1)" det(-J^iJi)) 

= (det Jip) (l + tr (tJ~lBi) + •••+ (-1) 2 " det (tJ^Bi)) , 

where tr (A) denotes the trace of a matrix A, and we used the well-known fact that 
the coefficients of the monic characteristic polynomial 

det(xl - A) = pa(x) - x n + a n _lir" -1 + ••• + a\x + ao, 
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are given by 

a n - k = (-l) k Y AjiA^ k= a 0 = (-l) n det(A), 

where {Aj }" =1 are the eigenvalues of A (see [ 8 ]). We define 

Ci, H {s) :=tr (Ji t o(s)Bi(s)) , C i>R (s,t) := t n p_j-i B .(l/t) - 1 - tC iyH (s), 
so that 


det = (det J i ,o(s))(l + fC' iiff (s) 


Consequently, we have 


(3.4) 


f fdH n 1 = f f(ai(s) + ti/(ai(s)))\detJ i)V (s,t)\ds 

JptXWiCT) JSi 

= f f(oti(s)+tv(ai(s)))\detJi t o(s)\\l + tCi iH (s)+ C itR (8,t)\ ds. 

JSi 


By Assumption 2.1, T is a compact C 2 hypersurface, and so the eigenvalues of 
Ji,o(s) and J-qB^s) are bounded uniformly in s e 5,;. Hence, for each 1 < i < TV, 
there exists a constant Ci such that 


|det - det Jj, 0 (s)| = |det J il 0 (s)| \tC i>H (s) + C R (s,t)\ < c t \t\. 

Moreover, by the compactness of T, there are only a finite number of Cj, and so we 
can deduce that there exists a constant c (that can be chosen independent of rj) such 
that, for all t e (-77,77), s e Si, 1 < * < N, 


max I det Ji n(s,t) - det Ji o(s)| < c\t \. 

l<i<N ’ 


Let {pi}Y be a partition of unity subordinate to the covering {Wi n r}(Y 1 of F. 
Consequently, by the diffeomorphism Q v , we observe that {/x^ o Pf 1 }^! is a partition 
of unity subordinate to the covering {pi{Wi nT)}^ of r t for t e (- 77 , 77 ). We define 

N N 

Ch(p) : = C R(P,t) : = Y (p) ^) ■: for per. 

2=1 2=1 

From the above discussion, we observe that C B (p) and C B (p,t ) are uniformly 
bounded in p eT and 


I tC H (p) + C R (p,t)\ < c\t\ for all \t\ < 77 . 

Then, from (3.3), for any / e L 1 (T) we have 

r N r N r 

/ fdU "- 1 =Y / pjd-H 11 - 1 = Y / |det J,;, 0 (s)| ds, 

Jr JWi<-\T JSi 
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and similarly from (3.4), for any / e L 1 (Tub I, (r)), 

f f(x) dx = f [ f dTL n ^ 1 dt = f V f (ni op~ l )fd'H n ^ 1 dt 

JTub"(r) Jr t J-r, i^JptiWinvy 

rv AL r 

= / Z/ Mi(a»(s))/(ai(s)+ti/(ai(s)))|det J ii o(s)||l + tCi,ij(s) + C'i,ii(s,£)| dsdt 

i=i “ / ‘S> 

= f Z f p i (p)f(p + tu(p))\l + tC H (p) + C R (p,t)\d'H n ~ 1 (p)dt 

J-r] JWinT 

= J J^f{p + tv(p))\l + tC H {p)+C R (p,t)\d'H n ~ 1 {p)dt. 


~v 

/- v Jr 
Hence, we can identify 


|det((Vp t ) T (Vp t ))| 2 (p,t) = |1 + tC H (p ) + C R (p,t )|. 

We summarise our findings of this section in the following: 

Lemma 3.1. Suppose Assumption 2.1 is satisfied. There exists c > 0 such that 
for any p < (c) -1 ; and any f e L 1 (Tub ?? (r)), we have 


(3.5) f /(x)dx=f [ f(p + tis(p))\l+tC H (p) + C R (p,t)\ d'H n 1 dt, 

Jiw>(r) J-r] Jr 

where 

(3.6) | tC H (p) + C R (p,t)\ < c\t\ for all \t\ < p. 

Consequently, 

T--r- [ |/(x)|dx< f f \f(p + tv{p))\dU n ~ 1 dt 

1 + cn J Tub^CD J-n Jr 


(3.7) 


/-rj 

1 


—— [ 

1 - cn Jt 


cp JTub 7 >(r) 


1 / 0*01 dx - 


3.2. Coordinates in a scaled tubular neighbourhood. In the subsequent 
convergence analysis, we will use a tubular neighbourhood whose width scales with 

, k 

e k for some 0 < k < 1, i.e., we consider X e := Tub c ?, (r). For this section, we take 
X s = Tub el? (r), i.e., k = 1, to derive some technical results. 

Let W c 1" be one of the open sets in a finite open cover of T, with associated 
local regular parameterisation a : S -* W n T and parameter domain S c R" _1 . We 
define the metric tensor Qo = ( 9 o,ij)i<%,j<n-i by 

(3.8) g 0 ,ij(p) : = d Si a(s) ■ d Sj a(s) for 1 <i, j < n- 1, if p = a(s) e T n W, 

with inverse metric tensor Qq 1 = ( 9 o)i<i,j<n-i■ Then, for a scalar function / : F -»■ R, 
the surface gradient of / on T at a point p = a(s) eT r\W for some s e S is defined as 
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One can show that the definition (3.9) does not depend on the parameterisation. 
There is also an alternate and equivalent definition of the surface gradient in terms 
of level sets, we refer the reader to [15, §2.2] for more details. 

Recall the signed distance function d defined in (1.1). We introduce the rescaled 
distance variable 

d 

z = —. 

£ 

For z e (— 77 , 77 ), we define a parallel hypersurface at distance ez away from T as 

(3.10) T EZ := {p + ezp(p) : p eT}. 

Let p(x) denote the closest point operator of x e X e as defined in (2.4), such that 

(3.11) x = p(x) + ezv(p(x)) for some z e (- 77 , 77 ). 

Then, by the injectivity of the closest point operator, we have 

(3.12) v(y) = u(p(y)) for y € T ez . 

For any scalar function / : X e -*■ R, we define its representation F e ( p , z ) in the 
(p, z) coordinate system by 

(3.13) F s (p,z) := f(p + ezv(p)) for p e T, z e (- 77 , 77 ). 

Locally, for x e X s n W, we have 

(3.14) x = G e (s,z ) := a(s) + Ezv(a(s)) for some s e S,z € (- 77 , 77 ). 

Then, we can define the local representation F e (s,z) of / in the (s,z) coordinate 
system by 

(3.15) F e (s,z) := F e (a(s),z ) = f(a(s) + szv(a(s))) for s e S,z e (- 77 , 77 ). 

Let (si,..., s n _i) € S and s n := z. Then, by (3.14), 

(3.16) d Si G e (s,z ) = d Si a(s) + Ezd Si v(a(s)) for 1 < i < n- 1, d Sn G e = £is(a(s)), 

and {d Si G e }™ =l is a basis of M" locally around V' sz . For 1 < i,j < n- 1, we define the 
metric tensor Q ez = (g e ,ij)i<i,j<n in these new coordinates as 

g e ,ij(p,z ) = (d Si a(s) + Ezd Si v(a(s))) ■ (d 8j a(s) + Ezd Sj v(a(s))), 

(3.17) g s ,in(p, z) = g e ,ni(p, Z ) = (d Si a(s) + £zd Si v(a(s))) ■ £v(a(s)) = 0, 

ge,nn(p,z) = (e!/(a(s))) • (£l7(a( S ))) = £ 2 , 

for p = a(s) e r n W, and we have used that d Si v ■ v = |i /| 2 = 0. 

Qez ~ (g l J)l<i,j<n denote the inverse metric tensor. Then, for any scalar 
function / : X £ -*■ R, with the representation F e (s , z) in the (s, z) coordinate system, 
we can express the surface gradient \7r sz Fs(s,z) on the parallel hypersurface r ez in 
local coordinates as 
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Similarly, we can define the surface gradient Vy'F e (s, z) on the original hypersurface 
T as 

(3.19) Vr F e (s,z)= Y Y 9 o ( a ( s )) d s i F e (s,z)\d Sj a(s). 

j =i \2=1 / 

Using (3.13) and (3.15) to switch to the (p,z) coordinate system, we can define 
the surface gradient X7re*F E (p,z) on the parallel hypersurface T ez via the formula 

(3.20) Vr ezF E (p, z) := Vr <**F E {oT l (p), z) for p = a(s) e T n W, 
and the surface gradient \/rF e (p, z) on the original hypersurface T as 

(3.21) Vr F E (p,z) := V r F E {a^ 1 (p),z) for p = a(s) e T n W. 

It is convenient to define the remainder both in local and global representations: 

(3 22) S/ S zF e (s,z) := Vr»F e ( 8 ,z)-VrF e (s,z), 

V ez F e (p , z) := Vr^F E (p, z) - Vr F e {p, z ), 

with the relation 


Ve Z F E (p,z) = \J ez F E {a 1 (p),z) for p = a(s) e T n W. 


We now state a result that decomposes the Euclidean gradient V/ into a surface 
component and a component in the normal direction: 

Lemma 3.2. Suppose that Assumption 2.1 is satisfied. Let f : X s ->• M be a C 1 
function with representation F e in the ( p , z) coordinate system and representation F e 
in the ( s,z ) coordinate system, as defined in (3.13) and (3.15), respectively. Then, for 
p € r, s £ S, z e (- 77 , 77 ) such that p = a(s), and x = a(s) + Ezv(a(s)) = p+ezv(p) € X s , 


(3.23) 


V/(x) = -u(p)d z F E (p,z) + V r ezF E (p,z) 

£ 

= -v(a(s))d z F E (s,z) + Vr^F e (s,z), 
£ 


where Vr«( - ) and Vr £z ( - ) are defined in (3.20) and (3.18), respectively. 

In addition, the remainders SJ EZ F E and S7 sz F e defined in (3.22) satisfy 

(3.24) \7 EZ F E (p,z) ■ v{p) = 0, S7 EZ F E (s,z) -v(a(s)) = 0, 

(3.25) S/ ez F e (p,z) = 0{e), S7 EZ F E (s,z) = O(e) as e 0. 


Proof. For the first assertion, we follow the proof given in [1, Appendix]. The 
equivalent result in two dimensions can be found in [23, Appendix B] and in [16]. 

Let / : X e -> i be a C 1 function. Recall the metric tensor Q EZ = (g e ,ij)i<i,j< n as 
defined in (3.17). If we denote the submatrix Q ez and its inverse Qf z by: 

Qez ~ (dJe,ij)l<i,j<n-l • Qez = (.9e ) • 


Then, we observe that 


Qez = 


/ Qez 0 \ £-1 _ / Q ez 0 \ 

l ° T r e2 i ° T / 
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where 0 e R" 1 is the zero column vector of length n — 1. Moreover, for f(x ) = 
F e (s(x),z(x )) 1 we have from (3.16) and (3.18), 

V/(z) = Y, [Y,9l J d Si F e d Sj G E = -d z F E d z G E + V r »F e = -d z F e v(a(s )) + Vr »F e . 

For the assertion regarding the remainder, let 

Cy (s) := d Si i/(a(s)) • d Sj a(s) + d Si a(s) ■ d Sj v(a(s)), 

’Dij(s) := d Si v(a{s)) ■ d Sj v{a{s)), 

and 


•- (1'/J ) 1 <i.j<7i- 1 ■ 7} — (H'ij ) 1 <i.j<ri -1 ■ 

so that by (3.17) and (3.8), 

5e,ij( a ( s )>^) = ffo,u(a(s)) + szCij(s) + (, £z) 2 T>ij(s ), 

Gez - Go + ezC + ( ez) 2 T>. 

A calculation involving the ansatz 

(3.26) 0 ” 1 = (Go + ezC + (^) 2 P ) _1 = G 0 1 + 5, 
will yield that 

£ = -{I + G 0 \ezC + (ez) 2/ D))~ 1 (Go l (ezC + (e*) 2 ®)^ 1 ), 

if Gez, Go and I+Go 1 ( ez C+(Ez) 2 V) are invertible. Here, I denotes the identity matrix. 

Since T is a compact C 2 hypersurface, all entries in the matrices Go, C, and T> are 
bounded. For a matrix H and e sufficiently small so that the absolute values of the 
eigenvalues of eH are less than 1, we have, 

(I +eHy l = I-eH+ e 2 H 2 

Hence, we can express 

GZl = Go 1 - ezG^CGo 1 + 0(e 2 ) as e - 0 . 

Consequently, looking at (3.18), and (3.19), we see that for s e S,z e (- 77 , 77 ), 

~ n—1 ___ 

\7 r e*F e (s,z) = Y 9 l J (a(s),z)d Si F e (s,z)d Sj G e (s,z) 
ij =1 

n —1 _ n —1 

. 00 _, = Y 9 o (u(s))d Si F e (s,z)d Sj a(s) + ez Y 9 o i^(s))d Si F E (s, z)d Sj v(a(s)) 

(3.27) i,j= 1 *,j=i 

n—1 _ 

-es ^ (^o 1 CC 7 o 1 )y(s) 5 Si F e (s,z) 5 S 3 a(s) + h.o.t. 

*d=i 

= V r F e (s, z) + V ez F e (s, z) = V r F e (s, z) + 0(e) as £ ^ 0, 
where h.o.t. denotes terms of higher order in e. By definition and by (3.12), 
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The analogous statements for the representation in the ( p , z) coordinate system 
can be obtain from (3.27) and (3.28) by applying the relations (3.20), (3.21), and 
(3.22). □ 

Remark 3.1. In two dimensions, for a local arclength parameterisation a = a(s), 
s e S c ]R, of T, we have that Qo = I, d s v - -nd s a, C = -2k, and Vr F e {s,z) = 
d s F e (s, z)d s a(s), where k is the curvature of T. Hence, we obtain from (3.23) and 
(3.27) that 


V/ = -d z F e v + d s F e d s a + Kezd s F £ d s a + <D(e 2 ), 
e 

which is consistent with [23, Appendix B[. 

Remark 3.2. For fixed e > 0 sufficiently small, such that dX e n <9f l = 0, we 
see that dX e has the same regularity as T. Using that C 1 (X e ) functions are dense 
in H 1 (X e ), we can extend the assertions of Lemma 3.2 to H 1 (X e ) functions, where 
(3.23) hold for weak derivatives, and d z (-) and Vr«(') are to be understood in the 
weak sense, and (3.24) and (3.25) hold almost everywhere. 

Using (3.5), (3.23), and a change of variables 1 >-*• ez, we have for any 0 < k < 1, 

f » , J /! 2 + l V /| 2 dx 

jTub £ ''(r) 

= k f r (\f\ 2 + \Vf\ 2 )(p + ti>(p))\l + tC H (p)+C R (p,t)\ dU 71 ' 1 dt, 

(3.29) £ _L r i 

= J ^ J r ~ \d z F e \ 2 (p, z)\l + szC H (p) + C R (p, ez) \ dOT~ x dz 

e 1-fc 

V 

+ XT X e(\F c \ 2 + |Vr-F e | 2 )(p, z) |1 + ezC H (p) + C R (p,ez)\ d H n ~ x dz. 

In addition, using that for e e (0,1], 

1 - ce k r] > 1 - cr/, 1 + c.e k rj <l + crj, 


we find that an analogous statement to (3.7) holds for Tub e 77 (T) for any 0 < k < 1: 


(3.30) 


k 

- f |/| dx < f f \f(p+ tv(p))\ d'H " -1 dt 

1 + dr/ JTub^»(r) IJI X £ ‘,ir 

- —~— f h |/| dx. 

1 - cn JTuh e ‘ , >(r) 


3.3. On functions extended constantly along the normal direction. Let 

/ e H l { T), and f e denote its constant extension off T in the normal direction to 
Tub r, (r), as defined in (2.6). A short calculation shows the following relation between 
V/ e and Vr/ (see [15, proof of Theorem 2.10] for more details): 


(3.31) Vf e (x) = (I - d(x)H(x))Vrf(p(x)) for x e Tub i 7 (T), 


where, H denotes the Hessian of the signed distance function d, and I denotes the 
identity tensor. Consequently 


(3.32) 


V/ e (y) = Vr f(y) for y e T. 
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Corollary 3.1. Let Assumption 2.1 be satisfied, and choose 77 > 0 with the 
diffeomorphism Q v : Tub ?, (T) -> Tx (- 11 ,rf) as defined in (2.3). Let 1 < q < 00 , 
f e W l ’ q (T), and let f e denote its constant extension in the normal direction offT to 
Tub f, (r), as defined in (2.6). Then, f e € lF 1 ’ ? (Tub ?7 (r)), and there exists a constant 
C >0, independent of f, such that 

II f L II-M(Tub 77 (r)) ^ Cl/IU^r), ||V/ c ||z /9 (Tub’’(r)) ^ Cl Vr/li^r)- 

Moreover, let f Ec denote the extension of f e from Tub ,; (r) to 12 by the method of 
reflection. Then f Ec € W l ' q (fii) and there exists a constant C > 0, independent of f, 
such that 

(3.33) \\f Ec \\Li(n)) ^ C1l/IU«(r)> II V/ Ec ||i,®(n)) ^ Cl Vr/IU^r)- 

Proof. Since r is a C 2 hypersurface, we see that 
(3-34) l-H’llc°( X ii b ''( r )) ^ IMIc 2 (Tub”(r)) < 00 • 

Let 1 < q < 00 , / e W 1,q (T) and let f e denote its constant extension in the normal 
direction off T to Tub T, (r), as defined in (2.6). Then, by (3.31), (3.5), (3.6), and (3.7), 
we have 

f \f e (x)\ q dx = f \f(p(x))\ q dx 
JTub'i(r) v ' dTub^cr) v v ' 

-Ur \f(p)\ q \1 + tC H (p) + C R (p,t)\ ATT 1 dt 
^C{cri)\\f\\ q Lq{TV 

f J v / e ( x )l ?dx =/' „. Urf(p(x))-d(x)H(x)V r f(p(x))\ q dx 
JTub’)(r) JTub’)(r) 

< C(1 + ||d||c 2 ( Tub ’i(r))) f \Vr f (p(x))\ q dx 

v v " JTub^r) 

< <7(1 + ||d||c 2 (Tub'<(r)))C'(c77)||Vr/||®, (r) - 
Thus, f e e W 1 ’ q (Tab ri (T)) with 

ll/ e ||iu 1 .<!(Tub’)(r)) ^ C(cp, |d||c 2 (Tub''(r)))||/||w 1 .<!(r)- 
For the case q = 00 , thanks to (2.6), we have 

II/ e II i~(TubV(r)) = ||/||L~(r)i 


and by (3.31), we have 

II V/ e ||i~(Tub 7 '(r)) $ C(crj, ||d|c 2 (Tub 7 '(r)))|| Vr/IU~(r)- 

By the extension theorem [18, Theorem 1, p. 254], we can extend f e to f Ec € 
with 

\\f Ec \\Li(n) ^ C'll/ e ||L‘!(Tub’’(r)) ^ C||/||L5(r)> 

II V/ Ec ||L9(n) ^ C|| V/ e ||z / q(Tub''(r)) ^ C|| Vr/||L®(r)! 
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where C is independent of /. □ 

The next lemma allows us to test with extensions of H 1 (T) functions as con¬ 
structed in Corollary 3.1 in the weak formulation of the diffuse domain approxima¬ 
tions. 

Lemma 3.3. Suppose Assumption 2.1 and 2.8 are satisfied. Let f e if 1 (T) and 
let f Ec denote its extension to f l as constructed in Corollary 3.1. Then, for all £ > 0, 

f Ec eH\Q,S e ), 


and there exists a constant C > 0, independent of f and e, such that 
\\f Ec \\o,6 s < C\\f\\ L 2 (r) , |V/ Bc ||oa < C'||V r /|| i2(r) . 


Furthermore, let b Ec and B Ec denote the extensions of the data b and B as mentioned 
in Assumption 2.6 and Remark 2.3. Then, as e -»■ 0, 


(3.35) 


[ s e b Ec \f Ec \ 2 dx [ b\f\ 2 dn n -\ 

Jn Jr 

[ 5 e B Ec Xlf Ec - V/ Bc dx - f BVrf • Vr/d'H" -1 . 
Jn ' Jr 


Consequently, as £ -»■ 0, 

(3.36) f 5 e \f E f dx - f |/| 2 d7T " 1 , f S e |v/ B f dx - f |V r /| 2 d U n ~ l 

J n Jr Jn 1 1 Jr 


Proof. We note that (3.36) follows from (3.35) if we consider b(p) = 1 with 
b Ec (x) = 1, and 5y(p) = Sij with b Ec (x ) = Sij, for p e T, x e Cl, 1 < i,j < n and Sij is 
the Kronecker delta. 

Choose r) > 0 and the diffeomorphism 0 17 as defined in (2.3). Thanks to (2.10) 
with <7 = 1 , 

(3.37) |<Mi°°(n\TW'(r)) 0 as e -*• 0, 

and so, for e € ( 0 , 1 ], 

(3.38) pe|i“(n\TW'(r)) ^ sup ||<5 e ||.L~(n\Tub T >(r)) = : Cs Up - 

£ 6 ( 0 , 1 ] 

Then, by (3.5) with a change of variables t >->■ es, (2.6), (3.6), (2.8), and (3.33) 
with q- 2 , 

J^e\f I - f Tuh n( r ^ e \f I dx + ll^e||L~(n\TW'(r))II/ llL 2 (fi\Tub’ ! (r)) 

< f ^ J r ^^j\f e (p + tu(p))\ 2 \l+tC H (p) + C R (p,t)\ dU n ~ x dt +C sup ||/ £;c ||^2 (n) 

77 

< f r $(s)\f e (p + esv(p))\ 2 \1 + £sC h (p) + C R (p,es)\ dU 11-1 ds + C sup \\f Ec \\ 2 L2(n) 

£ 

< f* J r S(s)(l + c£\s\)\f(p)\ 2 d H"- 1 ds + C sup ||/ Sc || 2 2 (n) 

< (1 + crj )|/|| 2 2 + C sup \\f Ec \\ 2 L 2^ < C\\f\\ 2 L 2( r y 
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Hence, f Ec e L 2 {VL,5 e ). 

Next, we have by the definition of b e , (2.6), (3.5) with a change of variables t >-» es, 
and (2.8): 

H 

f b e b e \f e f dx = f" f 6(s)b(p)\f(p)\ 2 \l + £sC H (p) + C R (p,£s)\ d'H n ~ 1 ds, 
J Tub r '(r) Jzi Jr 

[b\f\ 2 dU^ = [ [ 6(s)b(p)\f (p)f dH"-' ds. 

Jr Jr Jr 

Together with (3.6), and (2.9), we deduce that 


If 6 e b e \ff dx- [ b\f\ 2 AH™' 1 

L/Tub^rn Jr 


Xs(^,f) ^ dS H/Hi 2 (r) 

H 

Xa X^ S ^^l^^| 2 | £SC '' ff ^) + C, ' R ^’ eS )l d ^" _1 ds 

^ 2 <5(s)ds +ecC' 5>int J||6||i~ ( r)||/||i2 (r) ^0 as e ^ 0. 

e ’ e / / 


Then, the first statement of (3.35) follows from the above and 
f Seb Ec \f Ec \ dx < |< 5 e j| i ~(n sTub »( r ))!|fo' Ec | i =o(n)||/ £;c Ili 2 ( n ) -»■ 0 , as e -* 0 , 

where we have used (3.37). 

We point out that the proof for the assertions for the gradients follow similarly 
by using (3.31), i.e., 

l|V/ EC |lo,5 e ^ C'sup || V/ £7C llz,2(f2^ Xub T '( r >) 

H 

+ C(l^|c 2 (Tw>(r))) IIL <5(s)(l + ce|s|)|V r /(p )| 2 d'H n 1 ds 
< C sup ||V/ Bc || 2 2(n) + C(d,7?,c)||Vr/|| 2 L2(r) < C|| Vr/|| 2 2 (r) , 

which then implies that f Ec € H 1 ( 12,<5 e ). Similarly, using (3.5), (3.31) and the defini¬ 
tion of B e , and a change of variables t»£s, 

f 5 e B e (x)V f Ec (x) ■ S7 f Ec (x) dx 

JTub’J(r) 

= x ; x x ( j ) 6(p)[(jr “ Vr/(p)] ■ [(/ ■ tH)Vr/w] 

x|l + tC H (p) + C R (p,t)\ d H n ~ l dt 

?7 

= f: f r S ( s )B(p)[(I - esff)Vr/(p)] • [(/ - esff)Vr/(p)] 

x 1 1 + esCh ( p ) + C R (p, es)| d'H n ~ 1 ds 

77 

= : Xf f r s ( s ) B (p)Vrf(p) ■ Vrf(p)\l + esC H (p) + C R (p,es)\ d'H " -1 ds + 77, 
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where, by (3.34), (3.6), and (2.9), 

\R\ < £C(||d||c 2 ( Tub 7 ,( r )),c 77 ,C' < 5 j int)||#||L“(r)||Vr/||L 2 ( r ) -»■ 0 as e -*■ 0 . 
Thus, setting L := fJ 1 f r 5(s)B(p)Vrf(p)- Vr f(p) d l~L n ~ Y ds allow us to compute 


[ 6 e B Ec \7f Ec -Vf Ec dx-L + L- f BVrf-VrfdH 
J Tub^r) ' Jr 

f , v 3 /( s ) [ B(p)Vrf(p)- S7rf(p)d'H n ~ 1 ds 

V 

J-‘ S(s)B(p)Vrf(p) ■ Vr f(p) \esC H (p) + C R (p , es)| dH”" 1 ds 

2 d(s)ds +£cC 5 ,i n t| ||B||i~ (r )|| Vr/||i2 (r) + \U\ -> 0 as e -*■ 0. 
£ 5 £ ' / 


R\( 


Then, the second statement of (3.35) follows from the above and 


f 

J n 


nNTub”(r) 


6 E B Ec \7f Ec -Vf Ec dx < |<Mx,~ ( nx'iw- 


’(r)) 


I B 


Ecu 


Il°°( 0 ) 


llv/ 


Ec || 2 


Il 2 (si) 


0 , 


as £ -»■ 0 . □ 

Remark 3.3. For the double-obstacle regularisation, we use the fact that 5 e = 0 
on fl \ Tub^T) for rj > to deduce the same results. 

3.4. On the regularised indicator functions. Due to the boundedness of £ e , 
(2.7), and Lebesgue’s dominated convergence theorem, we have 

Lemma 3.4. Suppose that Assumption 2 .7 is satisfied, then for any g € L 1 (il), 

/ £,egdx -* g |n* dx, / (1 - dx 0 , / £ e gdx 0 as e -> 0 . 

J n* Jn* Jn* Jq\Q* 


Lemma 3.5. Suppose that Assumptions 2.1 and 2.8 are satisfied. Then, there 
exists a constant C > 0 such that, for all f e fU 1 - 1 (D) and all e e (0,1], 

^T|/|dx <(71/1^1.1(0). 

In particular, there exists a constant C > 0 such that, for all f e H X (D) and all 

£e(0,l], 

1/10,<5 e £ CI/llffRO)- 


Proof. Choose 77 > 0 and the diffeomorphism Q v as defined in (2.3). Let / e 
W 1 - 1 ^). Then, by (3.38), 

(3-39) / S e I/I dx < C SU p||/|| Li(f2\Tub T? (r)) ^ Csupll/llwi.Rn). 

7o\Tub r i(r) 

Using the diffeomorphism Q v , / iTubRr ) 6 W 1,l {T x (- 77 , 77 )). By absolute continuity 
on lines for W 1,1 functions (see [45, Theorem 2.1.4] or [33, Theorem 1, §1.1.3]) there 
exists a version of /, which we denote by the same symbol, such that for a.e. peT, 
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it is absolutely continuous as a function of t e (- 77 , 77 ). With absolute continuity with 
respect to t, we have for t € (- 77 , 77 ), 

f(p + tv(p)) = f(p)+ f -^ r /(p + 0'(p))dC- 

Jo dC 

Then, by (3.5), a change of coordinates t = £s, (2.8), and (3.7), 


f 4 |/| dx = f f -s(-)\f(p + tv(p))\\l+tC H (p) + C R (p, t)\ dH n 1 dt 

JTub T >(r) J-ri Jr e \eJ 


dC 


|l + iC'ir(p) + C , Ji(p,i)| dTT" 1 dt 


-XIXrM;)[ |/wl+ XlXr /foc) l 

^ (ll7o(/)|||i(r) + ^ |V/(p + CKp))I dC d'H”' 1 )^ <5(s)(l + &|s|)ds 

< (1 + £cC<5,int) (||7o(/)||L 1 (r) + ^(^ 77 ) || vy|| Z,i(Tub^( r )) ) 

< (1 + cCg, mt) (l|7o(/)||i!(r) + C(£ 77 )||/1 1 (Tub 1 '(r))) • 

Using the trace theorem [18, Theorem 1, p. 272] and (3.39), we have 
fje\f I dx < ((1 + cCs ,int )(Ctr + C(CT])) + 

Csup ) II/IIVU 1 - 1 (f2) 7 

where Ct r is the constant from the trace theorem. 

Note that the second assertion follows due to 


l|V(|/| 2 )||ii(n) ^ 2||/|| i 2(n)||V/||L2(fi) < 2j|/||^i (n) . 


□ 

Remark 3.4. For the double-obstacle regularisation, we directly obtain 

f 4|/| dx < (1 + cC s , i n t)(C' tr + C'(c 77 ))||/|| M /i,i(Q). 

-/Tub 2 " (r) 


Lemma 3.6. Suppose that Assumptions 2.1 and 2.8 are satisfied. For f e 
TU 1 , 1 (S 2 ), it holds that 

(3.40) [ S E f dx -*• f 7o(/)dH”" 1 as e - 0. 

Jn Jr 


Moreover, for rj > 0 sufficiently small, if f e 1 < q < 00 , or if f = C' 1 (f2) and 

q = 00 , then there exists a constant C > 0, independent of f and e, such that 


(3.41) 


f 4 /dx- fyoCOdTT - 1 

JTub^r) Jr 


< Ce 


1-7- 


l/ll 


W 1 ’i( fl)- 


Furthermore, if f e 2 < q < 00 with 70 (/) =0, or if f e C' 1 (fl) with f |r= 0 

and q = 00 , we have 


X 


Tub^cr) 


41/| 2 dx < Ce 2_ « l/l 


2 

vpb^n)- 


(3.42) 
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Proof. By the trace theorem and Lemma 3.5, the integrals in (3.40) are well- 
defined. Choose 0 < r] < (c) _1 and the diffeomorphism Q v as defined in (2.3), where c 
is the constant in (3.6). 

By (2.10) with g = 2, we see that 




£-<■0 £ Z 


Let C SU p j?? := sup ee ( 0jl i p-<5 (2) < oo. Then, for all e € (0,1], we have 


(?) 


< C F 2 
— ^sup,?) 0 ? 


and so, from the monotonicity of <5(|s|), the symmetry of S(s) about s - 0, and (2.9), 

ds = 2 L~ ■ ds s (iT ) 

^ 2^/Csi]p.//CVintc < Cs, 


(3.43) 


for some constant C > 0. 

Let q < oo and / € W 1,q {VL). Then, by the diffeomorphism Q 11 , absolute continuity 
on lines, and Holder’s inequality, we have 


(3.44) 


\f(p + £su(p))~ f(p)\ = \ f -^~f(p + £(v(p))d( 
I Jo dC 


= £ 


J q Vf{p + eC,v(p))-v(p)d( <e|s| q ^ |V/(p + £(v(p))\ q d( j’ . 


Then, by (3.44), (3.5), a change of variables t >-*■ es, Holder’s inequality, and (3.6), 
(3.45) 


I [ , SJdx- /' 7 o(/)dH n - 1 

< f R f r S(s)\x(^L 2 )(s)f(p + £sv(p))\l+£sC H (p) + C R (p,£s)\-f(p)\ d If' 1 ds 

< f R f r $(s)\f(p)\\l-X(^.2)(s)\l+£sC H (p)+C R (p,£s)\\ d W 1 - 1 ds 

+ f f <5(s)X(zn n.){s)\f{p + £sv{p))-f{p)\\l+£sC H (p)+C R {p 1 £s)\ dH” -1 ds 
^ ||7o(/)||i,i(r) X^) +£ c|s|) ds 

1 

+ £(l + crj) J^X(zrr)(s)S(s)\s\— ^ |V/(p + eC^(p))| 9 d£ dTT -1 ds 

< ||7o(/)||i 1 (r) (x R s(^, 2 )(s) +ec|s|) ds 

+ e(l + erf) |T|~JJ(s)(l + |s|)|j^ " |V/(p + eC^(p))| 9 dC dH" -1 j ds. 

|1 — ~ 

Here, we have also used that |s| r < 1 + |s| for all s e K and 1 < r < oo. By a change 
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(3.46) 



= £ « 


of variables e<C f, and (3.7), we observe that 

_1 

|V/(p + eC,u{p))\ q dC d7T _1 j 

1 

ui: \Vf(p + tv(p))\ q dC dH n <C(cg)£ «||V/|| L 1 ( Tub’J(r))) 

Together with (3.43), and (2.9), we see from (3.45) that 

I f $ef dx - f 7o(/)d"H fi 1 < eC'|7o(/)||_Li(r) + e 1 q C\ V/||L 9 ('iw>(r)) 
UTub^cr) Jr v ' \ \ n 

<Ce 1 9 |/||wi.i 




which is (3.41) for q < oo. 

For the case q = oo, let / e C 1 (fl). Then, by the fundamental theorem of calculus, 
we have 


(3.47) 


\f(p + esv(p)) ~f(p)\ =e\ f Mf(p + £C,v(p))-v(p)d(, 

I Jo 

<e|s| sup |V/(p+ £(v(p))\ < e|s| j| V/||c°(Tub r '(r)); 

Ce[0,s] 


and hence, by a similar calculation to (3.45) yields, 

I f 8 E f dx - f 7o (/)d'H n_1 

UTub^r) Jr 

< e C||7o(/)||i 1 (r) + £(1 + eg) |r| j| V/j|c°(Tub»(r)) M ds 

^ C ,e ll/llc 1 (f2)- 

Next, by (3.37), we see that 


[ 

Jn 


n\Tub r '(r) 


^ I/I dx < ||^ E ||i / ~(n N Tub , '(r))ll/llL 1 (n) ^ 0 as e -* 0. 


Together with (3.41) for q = oo, we have that (3.40) holds true for all C 1 (fl) functions. 
Let > 0 be arbitrary, then by the density of smooth functions, for any / € W 1,1 (J7), 
there exists g e C' 1 (fl) such that 

\\f ~ dWw 1 ’ 1 ^) < C- 

Then, by (3.40) for C' 1 (fl) functions, there exists £q(C) such that for e < £q(C)> 


I f S E g dx - [ "fo(g)dn h 
Jn Jr 


<C- 


By Lemma 3.5, and the trace theorem, we find that for e < £o(C); 
|^< 5 e /dx - J^'yo(f)dn n ^ 1 

< I f 8 e (f- g)dx - g)dH n ~ 1 

I Jn Jr 

^ C||/ ~ jllwr^n) + 1/ - jlU^r) + C 

< (C + C tr )\\f - j||rvn(n) + C ^ (C + Ctr + 1)C- 


+ f S e g dx - f 'yo(g)dK n 1 
Jn Jr 


(3.48) 
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By the arbitrariness of £, we conclude that (3.40) holds for any / e l4 Al ' 1 (fl). 

For assertion (3.42), we note that for / e C 1 (fl) with / |r= 0, a similar calculation 
to (3.44) yields for 2 < q < oo, 


\f(p + esv(p))\ 2 < e 2 I f |V/(p + £<>(p))I dC 
I Jo 

2 

<£ 2 |s| 2_i ^ \Vf(p + £(v(p))\ q dC Y ■ 

Then, a similar calculation to (3.45), and using (3.46) yield 
I [ 5 e \f\ 2 dx 

I JTub’»(r) 

TJ 2 

+ £sC H (p) + C R (p,£s)\6(s)\s\ 2 1 \V f(p +£(,v(p))\ q dC j dU n 1 ds 

2 

< £ 2 (1 + cp) |r| 1 ^ ^f r f Q ‘ f (p + £ Cv(p))\ q dC dTL"" 1 j (f R 5{s){l + \s\ + |s| 2 ) ds j 

<C£ 2 9 ||/||^/i,9(n)- 


Using the density of C 1 (S2) in W 1,q {VL), and arguing similarly as in (3.48), we see 
that assertion (3.42) holds true for any / e W 1,q {Tl), 2 < q < oo with 7o(/) = 0. 

While, for the case q = oo, let / e C 1 (fl) with / |r= 0. We have 


\f(p + esu(j>))- f(jp)\ = 


and so 


( 3 . 49 ) If <S E I/I 2 dx 

l-VTub^r) 


<e 2 ||/ll 


\f(p + £Sv(p ))\ 2 <e 2 |s | 2 IIV/ll 

C 1 (Tub 9 (r))(l + £r l) ^( s )l s l 

v' JR 


2 

C°(Tub^(r))5 


2 ds < Ce 2 l/l 


2 

CRfi)- 


□ 

Remark 3.5. The analogous assertions for the double-obstacle regularisation 
with p = £^ follows along a similar argument. For this case, we point out that (3.43) 
is not needed. 

We define 


#r,oW :={/€^ 1 (U): 7 o (/)=0}. 


This is a closed subspace of R 1 (U) under || • ||| ff i(o) and hence is a Hilbert space itself. 

Lemma 3.7. Suppose that Assumptions 2.1 and 2.8 are satisfied. There exists a 
constant C > 0 such that for all f e 0 (f2) and all £ £ (0,1], 

ll/llo,i<5 e ^ CII/IIhRO)- 


Moreover, 





ANALYSIS OF THE DIFFUSE DOMAIN APPROACH 


29 


Proof. Choose 0 < r] < (c) 1 an the diffeomorphism Q 71 as defined in (2.3). By 
(2.10) with q- 2, we have 


(3.51) 

Hence, for £ e (0,1], 


I g II L°°(r2\Tiib T? (r)) -*■ 0 as £ -* 0. 


I £ II L°°(f 2 \Tub T, (r)) — SU P II ~^e II L“(n\Tub”(r)) C'sup.-l ■ 

e<=(0,l] 

Let / e Up 0 (fi). By (3.42) with g = 2, we have 

f-5 E \f\ 2 dx< f -S E \f\ 2 dx + |^<M.L~ ( n\TW> (r)) ll/ll 
Jn e JTub^rr) e E v v " 


L 2 (r2\Tub 77 (r)) 


^ C||/llffi(a) + Csu P ,-i||/||l2(q) < C\\f\\ 2 H l(Q), 


where C is independent of / and e. Next, take g e C ,1 (f2) with g |r= 0. Then, from 
(3.49) and (3.51) we have 

In £^ e ^ 2 ~ C£\\g\\ci( n ) + ll;MXlT~(nxTub^(r))ll5llL 2 (n\Tub i '(r)) “*■ 0 as £ -»■ 0. 

Using a density argument similar to (3.48), we see that assertion (3.50) holds true for 
any /e^ 0 (9). □ 

As a consequence of (2.11), we have 


(3.52) 




and thus, for all £ e (0,1], and any / e L 2 {VL,f e ), we have 


(3.53) 


f |/| 2 dx < f L\f \ 2 

Jn £ Jn 


dx. 


We now introduce the following weighted Sobolev space: 

Definition 3.1. 

L 2 (T x R, S) := j/ : T x R ->■ R measurable s.t. XX S(z ) |/(p, z)\ 2 d H n 1 dz < oo|, 
with the inner product and induced norm: 

(f,g)L2(r* r,5) : = £ fs(z)f(p,z)g(p,z)d'H n ^ 1 dz, |/|l 2 ( r xR , 5 ) == (/,/)i 2 (rxRi5) , 
j m j r 

along with the identification, 

f = g o f(p, z ) = g(p, z) for a.e. p e T and a.e. z e {t e R : S(t) > 0}. 

In the next two lemma, we give compactness results for V e and H 1 ^,S e ). 

Lemma 3.8. Suppose that Assumptions 2.1, 2.1 and 2.8 are satisfied. Let 
{w e } e(E (o.i] c V e and {ru e } ee (o,i] c U 1 (H,^ e ) denote two bounded sequences, i.e., there 
exist constants C > 0, independent of s, such that 




+ \u 


lo,<s e 


<c, 
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Then, there exist u,w e H l {i 2*) such that u e |n* converges weakly to u in H 1 (PI*), 
w E |n* converges weakly to w in H 1 (Pl*), along subsequences, and for any ip e H 1 ( 12), 

(3.54) X 8 e u E ipdx. -* X 7 o(m) 7 o(^) d"H” 1 as e ->• 0. 

Proof. By Assumption 2.7, £ e > \ in Pi *, and so, 

(3.55) \\u £ |n* Iffun*') < 2 f £ E \u E \ 2 dx < 2 f £ E \u E \~ dx < 2 C for all e e (0, ll. 

v ' J n* Jn 

Hence, {u e |n*} e6 (o,i] c H 1 (H*) is a bounded sequence. Then, by the reflexive weak 
compactness theorem [18, Theorem 3, p. 639], there exists a function u € H x (12*) 
such that 


u e |n*-“- u in i7 1 (fl*) as e -*• 0, 

along a subsequence. The proof for the weak convergence of w e |n* to w in H 1 (12*) 
follows along the same lines. 

Choose 0 < r) < (c) _1 and the diffeomorphism Q v as defined in (2.3), where c is 

k 

the constant in (3.6). We consider the scaled tubular neighbourhood X e = Tub £ ^(T) 
for 0 < k < 1. Furthermore, choosing q = > 1 in (2.10), and using (2.8) and a 

rescaling, we obtain 

<5 e (x) < ^ 0 as 0, 

for all idl \ X e . Here, we have used that if x e 12 \ X e , then |d(x)| > £ k p. Thus, we 
deduce that 


(3.56) l<yL~ ( nsX*) -»■ 0 as £-»■ 0. 

Next, let ip € if 1 (12). Then, by the Cauchy-Schwarz inequality and the uniform 
boundedness of {u E } £€ (o,i] in V £ , we have 


(3.57) 


/ S e u e (p dx 
\Jn\xe 


^ ll ue llL 2 (nsA75 e )ll¥ , llL 2 (nN.Y75 e ) 

1 

<C'||4l2= (nxXe) ||^IU 2 (n) -*■ 0 as e 


0. 


Hence, it suffices to look at the integral over X e . Let U e (p,z) and <h £ (p, z) denote 
the representation of u E and ip in the ( p , z) coordinate system, respectively. For 
convenience, let us use the notation 


(3.58) 


r)c,k ■ gl-k ' 


Then, as u e is bounded uniformly in H 1 ^,^), by (3.52), (3.29), and (3.6), we have 




(3.59) 


f £ E |Vu e | 2 dx> f eCt5 e |Vu E | 2 dx 

JX* JXe 

f f £C i 8(z)^\d z U E \ 2 \l + ezC H (p) + C R (p,£z)\ dH 11 ” 1 dz 


C^l-cp) [ f x(- Ve , k ,n., k )(z)S(z)^\d z U E \ 2 (p,z)dn n 1 dz. 

J M J 1 c 
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Multiplying by e on both sides of the inequality allows us to deduce that 

(3.60) X(-r u,k,rie,k)( z ) d zU E (p, z) 0 in L 2 (T x R, S) as e 0. 

By (2.8), 5(0) > eS e (x ) for all icfl. As ip e i? 1 (fl), using (3.29), and (2.8) we have 

IMIffi(fi) ^ /, |V ^' dx ^ 1 

(3.61) 


f x -M 2 ^*mb eS '™ 2 dx 


X‘ 

1 - cq 


1 


)(z)6(z)-\d z <b e \ (j>,z)dH n dz, 


<5(0) JmJr X( ~ r ' e ’ fc:r ' E ' fc) - ' ' - e 

and thus, 

(3.62) X(-„ s , k , %lli )(z)3 z $ e (p^) 0 in i 2 (r x R,5) as e -* 0. 

Meanwhile, by Lemma 3.5, we see that 

IMIIu(fi) - ^IMIo,$ s - C f xe ^ e 1^1 dx 

>C(l-crj) f f x(.- Vetk ,ri e , k )(z)6(z)\^ e \ 2 (p,z)d / H n ~ 1 dz. 
~/jr o/r 

Similarly, the uniform boundedness of u £ in T 2 (f2,5 e ) gives 

C> ||u E ||o 6 > f S e \u E f dx 
JX £ 

* (!-«?) f fx(-ve,k,Ve,k)( z ) S ( z )\U e \ 2 (P,z)dn n - 1 dz. 

•J M J 1 


(3.63) 


(3.64) 


Hence, by the reflexive weak compactness theorem, there exists a function u e L 2 (F x 
R, S) such that 

(3.65) X(-r, e , k ,ri e , k )(z)U e (p,z) -u(p,z) in L 2 (T x R, 5) ase^O, 
along a subsequence. By (3.60) we can deduce that 

(3.66) d z u = 0 on J := {z 6 R : 5(z) > 0), 

and so u = u(p) in J. Indeed, for any 4/( p , z) that is smooth and compactly supported 
in J, we have from (3.60), 


(3.67) 


f f 6(z)(d z U e 'F)(p, z) dH n 1 dz -»■ 0 as e -»■ 0. 
J-ve,k Jr 


Then, by (2.9), and the smoothness of T, we have 


X 


r 

J(~Ve,k,rie,k)nJ Jr 5{z) 

^ C(cr), |r|) ||\I>|| L~ c (r x j) J 


\d>\ 2 (p,z)dU n - 1 dz 


Mf 

(-Ve,k,rie,k)nJ <5(z) 


dx < ■ 


S’(z) 

S(z) 


T( P: z)eL 2 (TxJ,S ), 


This implies that 
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and so by (3.65), as e -»■ 0, 

f f Xj(z)6(z)(u e (p,z)^Y^f{p,z)\dn^ 1 dz 
J-Ve,k J r \ 6{z) ] 

-»■ J J S'(z)(ud/)(p,z)dT-L n ~ 1 dz. 

Thus, we see that, as e -* 0, 

f Ve k [ 5(z)(d z U e y)(p,z)dH n - 1 dz 

■J-Ve.k Jr 

= -[' f Xj{z)5'{z){U e ^){p,z) + xj{z)d{z){d z ^U E ){p,z)dU n ~ 1 dz 

J -Ve,k 

~ fj Jr z ) + <5(z)(«c>;z\E , )(p, z) dH " -1 dz 

= Jj J^S(z)(d z u’d!)(p, z) dH^ 1 dz. 

But by (3.67), the left hand side also converge to zero as e -* 0. Hence, for arbitrary 
that is smooth and compactly supported in J, we have 

(3.68) J 8{z){d z u^){p, z) d"H ra_1 dz =0, 

which implies that d z u = 0 a.e. on J. 

To finish the proof, we will show that for all ip e H 1 (Cl), 

(3.69) J S e u e ip dx J uxo(p) d'H n ~ 1 as e -> 0, 
and then the identification 


(3.70) u - 70 (u) a.e. on T. 

First, we note that by (3.6) and the definition of <J> 6 (p, z) (see (3.13)), we have 
for a.e. (p, z) e T x R, 


(3.71) X(-r, e , k , Ve , k )(z)$e(p,z)\l + ezC H (p) + Cr(p,ez)\ 2 ^ 70 W)(p) as e -> 0. 
By Lemma 3.5, we have the uniform boundedness of the norm: 

f f 5(z) |$ e | 2 |1 +ezC H (p) + Cr(p,£z)\ d'H 71-1 dz 

(3.72) J ~ Ve - k Jr 

= f x Je \p\ 2 dx < C'||^||| r i (n) . 

By (3.40), we have 


f 5 e \p 2 \ dx -* f |7o(<^)| 2 d H" 1 as e -»■ 0. 
Jn 1 1 Jr 


Furthermore, from (3.56) we see that 
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Hence, we deduce that 


(3.73) [ 5 e \(p\ 2 dx -> f | 7 o((/p )| 2 d H n 1 as e -» 0. 

Jx e Jr 


JX 

In particular, we obtain the norm convergence: 

~Ve,k 


(3.74) 


[ f 5(z) |<h £ | 2 (p, z) |1 + ezCh(p) + Cr(p,£z)\ AT-C 1 1 dz 

J-V e,k Jr 

f S s \<pf dx -*■ f |7o(v?)| 2 d-H"- 1 = f f 5 (z)\j 0 (ip)\ 2 (p)d'H n ^ 1 dz, 
Jxe Jr Jr Jr 


as e -»■ 0 . 

Almost everywhere convergence (3.71) and uniform boundedness of the norm 
(3.72) imply weak convergence in L 2 (Y x R,<5) [7, Proposition 4.7.12, p. 282], To¬ 
gether with the norm convergence (3.74) yields strong convergence in L 2 (T x R 6) [7, 
Corollary 4.7.16 p. 285], I.e., as e ->■ 0, 


(3.75) 


X(-ve,k,Ve, k )( z )®e(P> z )\ 1 + £zC h(p) + C R (p,ez)\ z 
7o(t0(p) in i 2 (TxR,i5). 


Recall that Then, by (3.6), 

(3.76) sup \ezCh(p) + Cr(p,ez)\ < sup ce\z\ < ce k p, 

(p,z)erx(-r) e , k ,n eik ) (p,z)irx(- Ve:k ,ri etk ) 

and hence, using the identity 


yja- \fb = 


a-b 
yja + \/b 


we obtain, as e -*■ 0 , 


(3.77) 


ess sup 

(p,z)eTx(-n e , k ,Ve, k ) 


|1 + ezC H (p) + C R (p,ez) I 2 - 1 


C£ k l 7 


< C£ k p 0. 


\Jl + C£ k p + 1 
By the weak-strong product convergence, we have 

/ S e u E <pdx 
Jxe 

(3.78) = f [ 6(z)x(-r! e , k , T! e , k )(z)(U e <5>e)(p,z)\l+£zC H (p) + Cr(p,ez)\ d?T _1 dz 

v/T 

f K z ) [ u(p) r ro(p)(p)d'H n ~ 1 dz = f W7o(y>) dW 7 ‘ 

Jr Jr Jr 


jn—l 


as £ -»■ 0 . 


For the identification (3.70), we use that £ < 1 => r) e f, > p to see that the weak 
convergence of U e to u also holds in the restricted space L 2 (T x (— 77 , 77 ), <5), which is 
equivalent to restricting to Tub er, (T). Indeed, for 'k € L 2 (T x R, S), we have 


(3.79) 


Jj(z)(U e -u)d> d'H n ~ 1 dz 

f fs( Z )(U e -u)(x(^r,)(zm dTT - 1 dz 

jM j1 
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as e -* 0. Let 



Then, by (3.6), we find that 

XIX b{z)\lo(v)(p)\ 2 \\ l +£zC H (p) + C R (p,£z)\-l\ d n n 1 dz 

X T] 

<5(0 M dz < C’(c,?7,C'5 i i n t)||7o(^)|li2 (r) e^ 0 as e ^ 0, 

and thus we deduce the following strong convergence: 

^3 8Q ) X(-r,,n){z)lo ( 0(0 |1 + ezC H (p) + C R (p,ez) | 5 

-»■ 7o (v)(p)X(-ri, v )(z) in L 2 (r xR,d), as e 0. 

By the Cauchy-Schwarz inequality, (3.6), (2.9) and Lemma 3.5, 

f f r ^(z)\'y 0 (u)(p)'y 0 (ip)(p)\\l + ezCH(p) + C R (p,ez)\ cl"H”" 1 dz 

- x; x <5(0 l7o(0(0l l7o(OCOI (1 + ce|z|) d'H n 1 dz 

< (1 + c£COnt)|7o(OIU 2 (r)|7o(OIU 2 (r) < oo. 


Recalling that, by definition (see (3.13)), <h £ (p, 0) = 7o(0; we can now compute 
(writing Ch = C R (p) and C R = C R (p,ez) in the following): 


|J f r 6(z)((U e $ e )(p,z)-(lo(u)jo( l p))(p))\l + £zCH + C R \dH n 1 dz 

< \f f r 6(z)(U £ (p,0)-~/ 0 (u)(p))"/o(v)(p)\l+£zC H + C R \ dU 71 ' 1 dz 
(3.81) + |/’'j/(0(^~ -^(U £ $ e )( P ,() d^ll + ezC H + C R \ dH™- 1 dz 

< (1 + cCi^nt) | 7 o(m £ ) - 7o(fi)| l7o(OI d'H n ~ l 


(i +-°i)f_yp)f r f_ 


—(U : ‘l/)(]>. C) 

dC 


dC d H n ^ dz. 


The first term on the right hand side converges to zero by the strong convergence 
of 70 (u e ) to 7 o(u) in L 2 (T) (see Theorem B.3). For the second term, using the 
monotonicity of 5(|s|), and the definition of C v , 


X>>XrX 


"HO 
V < 5(0 


dC 


(C/ £ 4> e )(p, 0 


dC d'H n ~ 1 dz 


(3.82) < f r £ S(0(\d c U £ ( P , OMp> 01 + \U e (p , 0 d c * e (j>, Ol) dC d U n ~ l 


< 5(0 

< C\d z U E ||i2( r x(_ T7 , T7 ) J 5)II||i2( r x(- 77 ,7 7 ),5> 
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as e -*■ 0, where we used (3.60), (3.62), (3.63), and (3.64) to show that 

\\d z U e \\ L 2 (rx(-i7,i7),5)ll < ^E||L 2 (rx(-r ? ,r ? ),5) 

+ l^ £ ||L 2 (rx(-^,^),5)||9. < i) e || L 2 (Tx(-ri,r)) ,5) 

^ C|| ( Plu 1 (f2)IX(-'? e ,fc,^.fc)( Z )^^ E |L 2 (rxK,<5) 

+ C , ||u e ||o,5 e |X(-r ;e , fc ,r ;e , )! )( Z )^ $ ellL 2 (rxR,5) -* 0, 

as e -* 0. As a consequence of (3.81) and (3.82), we have 

f fs(z)((U e ^> e )(p,z)-('yo(u)jo(tf))(p))\l + £zC H + C R \d'H n ^ 1 dz 
(3.83) J-v -'F 

-»■ 0 as e -»■ 0. 


But, by (3.6) and (3.80), 

f J r d(z) 7 o(u)(p)'yo(tp)(p) |1 + ezC H + C R \ d 'H n ~ 1 dz 

->■ J J^<5(2)7o(u)(p)7o(^)(p) dU^ 1 dz = C v J^'Yoiuhoiv) dH 71 ' 1 as £ -* 0 . 

Together with the weak convergence of U e to u in L 2 (T x (— 77 , 77 ), <5) (see (3.79)), we 
have 

f f 5(z)((U e $ e )(p, z) - (jo(u)"/o(p))(p))\l + £zCh + Cr\ d'H n_1 dz 

J -77 Jr 

^ X,, f r S ( z )(u(p)-lo(u)(p))io(<p)(p)dH n - 1 dz = Cn J^u - jo(u))-/ 0 (p) dW 1 ' 1 , 

as e -»■ 0. Hence, by (3.83), we have 

0 = C v J (u - 7 o(u)) 7 o(<^) dH " -1 for all p e i7 1 (H), 

which implies that u = 70 ( 77 ) a.e. on T. □ 

Remark 3.6. For the identification (3.70), we restricted to the tubular neigh¬ 
bourhood Tub 6 ,, (T), so that in the ( p,z ) coordinate system, we have z e (- 77 , 77 ). This 
is needed in (3.82) where we estimate the fraction from above. Otherwise, we 
cannot deduce that the right hand side of (3.81) converges to zero as £ -* 0. 

Remark 3.7. For the double-obstacle regularisation, we choose X e = Tub £ ^(T) 
and we will have U e converging weakly to u(p) in L 2 (Tx (~?,?)>^)- The estimate 
(3.76) becomes 

7T 

sup I £zCh{p) + Cr(p,£z)\<C.£-. 
per,« 6 (-f,f) 2 


A similar argument with the above elements will show (3.69), and we restrict to 
Tub ET (T) in order to show (3.70). 

Corollary 3.2. Suppose that Assumptions 2.1, 2.1 and 2.8 are satisfied. Let 
{w e } E€ (o.i] c bV e be a bounded sequence, i.e., there exist a constant C > 0, independent 
of £, such that 


lUe 


+ fur 






IKIIoa ^ c - 
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Then, there exists w e Hq(Q*) such that w e |n* converges weakly to w in 

Proof. For e € (0,1], we see that ||i(J e ||§ < Ce < C, and thus, by Lemma 3.8, 

there exists a function w e H 1 {PL*) such that, for all ip e H 1 (Q), 

w e |n* w in iL 1 (fl*) as e -* 0 , 

(3.84) J~ S e w e ipdx -*■ X 7 o(w) 7 o(^)d'H n 1 as e ->• 0 , 


along a subsequence. 

On the other hand, we have by Lemma 3.5, Cauchy-Schwarz inequality and the 
uniform boundedness of {'w e } e€ (o,i] in W e : 



f 'z8 e w e ipdx 

\Jn \fe 


I f dx ) 2 ( f eS e 

\J n e ) \J n 



£ £ Cy/e\ip\ H i(Q,-) -* 0 as e -* 0 . 


Upon comparing with (3.84), we see that 


d TL n 1 = 0 for all ip e i? 1 (0), 


and so, 70 (ui) = 0 a.e. on F. □ 

Lemma 3.9. Suppose that T is a C 3 compact hypersurface and that Assumptions 
2.1 and 2.8 are satisfied. Let {t’ £ } e€ (o,i] c F^ 1 (U,<5 e ) be a bounded sequence, i.e., there 
exists a constant C > 0, independent of e, such that, ||u £ ||^ $ < C. For 1 < i,j < n, 
let g Ec £ L 2 (Ll), b Ec e and b Ec e L°°(fi) denote the extensions of the data 

g e L 2 (T), be L°°(T ) and bij e L°°(T ) as mentioned in Assumption 2.6 and Remark 
2.3. 

Then, there exists v e H 3 ( T) such that, 


(3.85) 



and for any ip e H 1 (f2) and if) e /f 1 (r) with 1 p Ec e H 1 ( fl) as constructed in Corollary 
3.1, 


(3.86) 


f S e b Ec v s ipdx -* f bv^o(ip) dTL n 1 , 

J n Jr 

f 5 e B Ec X/v e • S/ip Ec dx - f SV r U • VrV’ d'H n ~ 1 , 
Jn Jr 


as e -»■ 0. 

Proof. Choose 0 < p < (c ) -1 and the diffeomorphism O ’ 7 as defined in (2.3). Let 
X s = Tub e ’'(r) for 0 < k < 1. By (3.56), and the uniform boundedness of {i> £ } e e(o,i] 
in H 1 (Cl,S e ), 


\[ 

Jn 


Q.\X e 


eV g 


Ec 


dx 


* ||« e ||o,«. II41l- ( nv!t.)||0 C iU 2 ( 0 ) 0 as er 0 



0 as e -»■ 0 , 


and 
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for any p e if 1 (51). Similarly, by (3.56), and the uniform boundedness of {u e } e6 (o,i] 
in H 1 { 51,<5 e ), for any ip Ec e H 1 { 51) constructed from ip 6 77 1 (T) as in Corollary 3.1, 


I f 

I Jsu 


/£2\A e 

< w Ec 


6 e 13 Ec S7v e -Vip Ec dx 


(n) IIV| o.5 e I S e 1 ^ xe) | Vip L || L2( n > 


0 as e 0. 


Hence, it is sufficient to restrict our attention to X e . Then, invoking the (p, z) coor¬ 
dinate system, and using (3.29), and (3.30), 


(3.87) 


C'>||u £ ||?, 5e > f <UH 2 + |VT £ | 2 )dx 

e J A' £ 

> (1 -crj) S(z) (\V e \ 2 + 1 \d z V e \ 2 + |Vr« V B | 2 ) (p, z) d U n ~ x dz. 


Hence, there exists a function v 6 i 2 (r x R, 5) and a vector-valued function Q e 
(T 2 (T x R, 5))™ such that, as e -*■ 0, 

X(- Ve , k , Ve , k )(z)d z V £ (p,z) 0 in L 2 (T xK,(S), 
X(- Ve , k , Vs , k )(z)V s (p,z) -v(p,z) in L 2 (T x R,5), 
X(- Ve , k ,ri e , k )(z)Vre*V E (p,z) - Q(p,z) in (L 2 (r x R, 5))”. 

Moreover, we can deduce, via a similar argument to the derivation of (3.68), that 
d z v = 0 and v = v(p) a.e. on J. 

We claim that Q = Vru(p). As T is a C 3 compact hypersurface, and so for a 
finite open cover {Wi nT}fi 1 ,Wi c R” of T, we have local regular parameterisations 
cii: Si -»■ Wi nT. Let {pijY denote a partition of unity subordinate to {W,}^. Take 
Y e C“(T x R), and for any 1 < r < n, let Q r (s,z), V e (s,z), T(s,z), VreT £ (s,z) 
denote the representation of Q r , V e , Y, Vr«k E in the (s, z ) coordinate system. Then, 

f f 6(z)(Q r Y)(p,z)d'H n ~ 1 dz 

Jr Jr 

N r r _ 

= Z/ / M*(ai(s))(5(2;)(Q r y)(s,z) |det J i>0 (s)| ds dz. 

Jr JS; 


Moreover, since T is C 3 , the normal i/ and the components of the inverse of 
the metric tensor Q o on T both belong to the class C 2 . From (3.17), we infer that the 
components g l e J of the inverse of the metric tensor Q ez are C 1 functions. So, for any 
1 < r < n, we have by (3.14) and integration by parts, 

(3.88) 

f pi(ai(s))5(z)(VrezV e ) r (s,z)Y(s,z) |det J i>0 (s)| ds 

JSi 

n—1 ~ 

IM(ai(s))S(z) Y 9 J s r (s, z)d Sj V e (s, z)d Sr G e (s, z)Y(s, z) |det Jj, 0 (s)| ds 
3=1 

n—1 _ _ 

S(z) Y 9 Sj (gi r (s, z)/Ui(ai(s))Y(s, z) |det Ji,o(s)| d Sr G e (s, z)) F £ (s, z) ds 
3=1 

Note that, by the regularity of T and Y, 

d Sj (gi r [HY\detJ ifi \d Sr G s ) e L 2 (S t x R ,8). 
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Hence, integrating with respect to z, summing from i = 1 to N and passing to the 
limit as e -*• 0 in (3.88), we obtain 

Y f f IM(cti(s))5(z)((Vre*V e ) r Y)(s,z) |det J i| 0 (s)| ds dz 

i= 1 J-Ve,k JSi 

A ^ r r n —1 _ 

-»■ ~Y / / s ( z ) Y {g : o r ( s )^i( a i( s ))^( s ^ z ) |detJ i| o(s)| 0 Sr a i )tJ(s)ds dz 

i=l “ /R j= 1 

iV r r n—1 __ 

= Z / / Aij(«i(s))<^(^) Z So r (s) 0 g i u(s)c? Sr a! i (s)Y'(s, 2 ) l det J r i,o(s)| ds dz 

i=i JR J=1 

= f f 5(z)(Vrv)r(p)Y(p,z)d'H n ~ 1 dz. 

Jr Jr 

Meanwhile, by the weak convergence of X(-r ls k ,ri e k )( z )(^r^V e ) r to Q r , we have 


SL rVe,k r _ 

Y / /r l (a l (s))((Vr«’l /e ) J .y)(s, 2 :)]det Ji,o(s)\ ds dz 

i=l J Ve,k JSi 

N r r _ __ 

V / / fj,i(ai(s))S(z)(Q r Y)(s,z) |det J lj 0 (s)| ds dz as £ -»■ 0. 
Jr Js 4 


Transforming to the (p, coordinate system and equating the limits leads to 

f f S(z)(Q r Y)(p,z)d'H n ~ 1 dz = f f (V r v) r (p)Y(p,z) dJT ’ 1 dz, 

and since Y is arbitrary, we have that Q r = (VrtJ) r . 

Let g Bc , b Ec and B Ec denote the extensions of the data g , b and B as mentioned 
in Assumption 2.6 and Remark 2.3. We note that by the definition of the extensions: 

g E \x) = g\x) = g(p(x)), b Ec (x ) = }f(x) = b(p(x)), B Ec (x) = B e (x) = B(p(x)) : 

for any x e X e . To show (3.85) it suffices to show the strong convergence in L 2 (T x 

K,<5): 

(3.89) X(.-r, e , k ,r, Stk )(z)g(p)\l + £zC H {p) + C R {p,£z)\^ -> g(p) as e -* 0. 

We use the method in the proof of Lemma 3.8. Note that (3.89) holds pointwise for 
a.e. p e T and zel. Moreover, by Corollary 3.1, we have uniform boundedness of the 
norm: 

f xe Ss\g e | 2 dx < C'||p|| L 2 (r) . 


Furthermore, we may appeal to (3.36)i to show norm convergence, since the proof 
also works for the extension of a function in L 2 (T). Then, arguing as in the proof 
of Lemma 3.8 yields the required strong convergence. Using the weak convergence 
X(-n e k ,ve k)( z )V e to v i n l2 ( t x IK.,<5), we obtain by the weak-strong product conver¬ 
gence: 


f 5ev e g e dx 
Jx e 

[ f $(z)X(- Veik ,ri e , k )( z )V e (p,z)g(p)\l + ezCh(p) + Cr(p,ez)\ dU n _1 dz 

jk jl 

f f 5(z)v(p)g(p)d'H n ~ 1 dz = [ vgdH n_1 as e -> 0. 

Jr Jr Jr 
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We note that the proof of (3.86)i is similar to that of (3.69). More precisely, we 
replace u E with v e , and u with v. Using that b Ec (x) = b(p(x )) for all x € X s , and 
(3.77), we can deduce the strong convergence 

b Ec (p + ezv(p)) |1 + ezC H (p) + C R (p, £z)\ 1 
= b(p) |1 + £zC h {p) + Cr(p,£z)\* -* b(p) in L°°(r x R,<5). 


Together with (3.75) and the weak convergence of X(-r) e k , Ve k ){ z )V e to v in L 2 {V x 
R, 5), we obtain by the weak-strong product convergence: 


/ 5 e b Ec v e ipdx 
Jx £ 

= f [ $( z )X(-ri e , k , Ve , k )( z )(V e $e)(p,z)b(p)\l + £zC h (p) +Cr(p,£z)\ dU n ~ l dz 
j m -J r 

f f S(z)v(p)'yo(<p)(p)b(p) d"H" _1 dz = f bvxo{p) d'H n ~ 1 as e -»■ 0. 

Jr Jr Jr 


Next, recall that for a given ip e JJ 1 (T) with extension ip Ec e H 1 { S2) as constructed 
in Corollary 3.1, within the scaled tubular neighbourhood X E , ip Ec is simply the 
constant extension ip e in the normal direction off T, as defined in (2.6). The same 
is true for the extension of the data B , i.e., B Ec (x) = B(p(x)) for all x £ X E . By 
Corollary 3.1, we have that ip e £ H 1 (X E ) and by (3.31) we have X7ip e (x ) • i/(p(x)) = 0. 
Moreover, let 4U denote the representation of ip in the (p , z) coordinate system. Then, 
by (3.23) and (3.24), we see that, for x = p + £zu(p) £ X E , 

(3.90) 0 = S7ip e (x) ■ v(j>(x)) = -<TT £ (p, z) + V r « ^ E (P, z ) • v(j>) = -d z ^ e {p, z), 

£ £ 

i.e. dz'kc - 0. Together with (3.24), we compute that 


(3.91) 


f 5^B e XJv e ■ Vip e dx 

[' hk f S(z)Vr^V E ■B(p) T Vr^e\l + £zC H (p) + CR(p,£z)\dH n - 1 dz, 
J —rip t Jr 


where B T denotes the transpose of B. 

With the weak convergence X(-r) £ k ,ve k )( z )Vr «U e (p, z ) to Vr«(p) in L 2 (TxR, 5), 
in order to show 



S E B e W E ■ Vip e dx 


J^BVrv-VripdH 71 - 1 


as £ -*■ 0 , 


it is sufficient to show the following strong convergence result: As e -*■ 0, 

^3 92 ^ X(-r, e , k:Ve , k )(z)B T (p)Vre^ E (p,z)\l+£zC H (p)+C R (p,£z)\ 3 

B T (p) Vr^(p) in L 2 (T x R, 5). 


By (3.13), (3.6), and (3.25), we have for almost every p e T, 2 e R, as e -* 0, 
X(-rt e , k , Ve , k ){z)B T {p) Vr«^ e (p, z) |1 + £zC h (p ) + C R (p,ez) I 3 

>3 T (p)Vrtp(p)- 


(3.93) 
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By (3.90), (3.23) and a slight modification to the proof of Corollary 3.1, we have 
uniform boundedness of the norm: 


(3.94) 


k) (z) |K T V r e--4' £ | |1 + ezC H (p) + C R (p,£z)\ dU n 1 dz 


lli 2 (r)- 


LL s(z)x <-'" 

5 e \(B T ) e vr\ 2 dx < C||B|| i - (r) ||VrVl|] 

Lastly, by a slight modification of the proof of (3.35)2 and using that 
f n xe S e \(B Ec ) T vy Ec \ 2 dx < !|4||L~(asA-)||B Bc |U~(n)||V^ c ||i 2(n) - 0 as e - 0, 
we can deduce 


[ 5 e \B T 'V'ip e \ dx -*■ f |6 T VrV , l d"H" 1 as e 
Jx * 1 1 Jr 1 1 


0 , 


and hence we obtain the norm convergence: As e -» 0, 


/ Ve,k r ,2 , 

/ S(z) \B T (p)Vr^^e(p,z)\ \l + £zC H (p) + C R (p,£z)\d'H n dz 
Ve,k Jr 

(3.95) = [ S e \B T vy e \ 2 dx -» flMdr 1 

Jx £ 1 1 Jr 1 1 

= f f 5(z)\B t (p)(V rip)(p)\ 2 d'H n ~ 1 dz. 

Jr Jr 


Then, arguing as in the proof of (3.75), we use (3.93), (3.94) and (3.95) to deduce the 
required strong convergence (3.92). □ 

Corollary 3.3. Suppose that Assumptions 2.1 and 2.8 are satisfied. Let g Ec € 
L 2 (Ll) denote the extensions of the data g e L 2 (V) as mentioned in Assumption 2.6 
and Remark 2.3. Let e H 1 ^). Then, as £ -»■ 0, 

(3.96) f S e g Ec ipdx f gj 0 (tp) d H n ~ 1 . 

Jo Jr 


Proof. Choose 0 < p < (c) 1 and the diffeomorphism Q v as defined in (2.3). Let 
X £ = Tub^T) for 0 < k < 1. By (3.56), 

\fn x J e ^9 EC ^ ~ IIML~(n\xe)lls £c ||L 2 (n)MI.L 2 (n) -*■ 0 as e -*■ 0. 


Hence, it is sufficient to restrict our attention to X s . Invoking the (p, z) coordinate 
system, we note that by definition g Ec (x) - g e (x) = g(p(x)) for all x € X e . Moreover, 
(3.71), (3.72), (3.74), and (3.75) are still valid. Hence, by (3.75) and (3.77), we have 

f ^eJV^dx = f ’ [ 6(z)g(p)$ £ (jp,z) |1 + £zC H {p) + C R (p,£z)\ dH n ~ l dz 
Jx* J-Ve,k Jr 

f f b(z)g(pYi Q {p){p) dTL n ~ l dz = f g ~/ 0 ( v ) dU n ~ l as e -*• 0. 

Jr Jr Jr 


□ 


4. Proof of the main results. 
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4.1. Proof of Lemma 2.1. The weak formulation of (RSIH) is: Find w e 
H 1 (12*) such that for all ip e H 1 ( fl*), 


(4.1) 


/ A\7w-\7(p + awipdx + / P^Aw^^Ap) d'H n 1 
Jn* Jr 

/ fip - AVg • V<p - agpdx. 

Jn* 


By the Lax-Milgram theorem, for each fj > 0, there exists a unique weak solution 
vjP e H 1 ( f2*) to (RSIH). Moreover, w^ satisfies the following estimate, with the 
constant CA.a independent of /3: 


min( 6 » 0 , 6 » 2 )||u; / 3 ||^i (n * ) + 2 / 3 || 7 o(w /:, )||i 2 (r) < 


/3U|2 


min ( 0 0 , 6 2 ) 


From this, we observe that {w^}p is bounded uniformly in H Hence there 
exists a subsequence (pi)i<M -* 00 such that 

W A ^ w ° e as Pi -* 00 . 

Moreover, from the a priori estimates on {w^}p, we see that 

|| / /3\|| ^ ^(^0; @2 1 /*, A 1 U, £/) 

ll7o(wC)llz, 2 (r) <- — -» 0 as p-> 00 . 

Compactness of the trace operator (Theorem B.3) implies that 70 (u>^) converges 
strongly to 70 ( 10 °), and subsequently the norms also converge. Thus, 

0 = lim |7o(w / 3 i )iU 2 (r) = |7o(^ 0 )IU 2 (r)i 

Pi—> 00 

and so 70 ( 10 °) = 0 a.e. on T. Since T is C 1 , by the characterisation of zero trace 
Sobolev functions [18, Theorem 2, p. 259] we have w° e Hq(CI*). Testing with 


<p e Hq(Q*) in (4.1) and passing to the limit P ■ 


we see that w° satisfies 


/ AV w° ■ V p + aw° ip dx = / ftp-AVg-Vp-agp dx, 

Jn* Jn* 

for all <p e Hq(Q*). Hence converges (along a subsequence) to the solution of 
(DSIH) weakly in H 1 (f2*) as p -*■ 00 . Since (DSIH) is well-posed by the Lax-Milgram 
theorem, the aforementioned convergence applies to the whole sequence. 

4.2. Well-posedness for (CSI). We consider the product Hilbert space and 
associated inner product 

X := x H\T), 

((ui,v 1 ),(u 2 ,v 2 ))x := L u\u 2 + Viti • Vu 2 dx + J~ iqi »2 + Vrih • Vri^ d VA 1 . 


Let 


acsi((u,v),(ip,i/j)) := a B (u,ip)+a s (v,ip) + Kl s (v - 70(11), i/> - 7o(<0)), 
lcsi((g>,ip)) ■= lB(f,y>) + Pls(g,i>), 
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where a^( as( ^s(v) and Zg(•,•) are as defined in (2.1) and (2.2). The weak 
formulation for (CSI) is: Find (u,u) e T such that for all e ff, 

acsi((u,v),(<p,ip)) =l C si((‘P,i’))- 
By the root mean square inequality, i.e., 

(4.2) a+b< V2,Va 2 + b 2 , 
one can show that 

Vcsi{{vA’))\ ^ l/|z, 2 (n*)IMlL 2 (n») + 1^1 z. 2 (r)Ill'lli 2 (r) 

(4.3) < (|/|U 2 (n*) + ||s , ||i. 2 (r))(IMU 2 (n*) + ||V , IU 2 (r)) 

^ v^(||/||i 2 (n*) + IIS'IIz, 2 (r))II(V?,■»/>)II at, 

and 

Jj(u-7o(u))(V’-7o(v?))l d U n ~ x 

< \\v - 7 o(w)||i 2 (r) \\ip - 7o( c ^ 7 ) II z, 2 (r) 

^ (IMU 2 (r) + Ctr|w||#i(n*))(||V , ||L 2 (r) + CtrlMIfl-^n*)) 

^ + C?r)K u ’V)\\xKv,ip)\\x, 

and so 

\acsi((u,v),(<p,ip))\ ^ (CA,B,a,b + 2K(1 + Cl))\\(u,v)\\x\\((p,^)\\x- 
Moreover, 

acsi((u,v),(u,v)) > min( 6 »o, 6 > 2 )||u|^i (n * ) +min( 6 »i, 6 ' 3 )|v|^i (r) + fs:||u- 7 o(M)||| 2 (r) 
> (min 9i)\\(u,v)\\%. 

I 

By the Lax-Milgram theorem, there exists a unique weak solution (u,v) € X to (CSI) 
such that 

II(m,u)|U < C'( 0 i)(||/||i, 2 (n*) + ||fl , IU 2 (r))- 

4.3. Well-posedness for (CDD). We consider the product Hilbert space and 
associated inner product 

X s :=V s xH\n,6 e ), 

{(u 1 ,v 1 ),(u 2 ,v 2 ))x e := J (€e + 5 e )uiU2 + €eVui • Vrt 2 + 5 e vit > 2 + 5 e Vvi • Vw 2 dx. 

Let 

(4.4) a C DD((u,v),((p,ijj)) := a e B (u, tp) + a%(v, tp) + Kl e s (y - u, tp - tp), 

(4.5) Icdd((<P,1>)) := l%U Ea ,<p) + Pls(9 Ec ,^), 

where a s(' 7 ')) a s(v)i (•,•) and Z|(-,-) are as defined in (2.15) and (2.16). The weak 
formulation for (CDD) is: Find (u e ,v E ) e X e such that for all (p,tp) e X s , 

a C DD((u E ,v s ),(^)) = 1 C dd((<P,^))- 
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A similar calculation to (4.3) involving the root mean square inequality (4.2) applied 
to 1 C dd yields 


\Icdd((v,tI>))\ <V2(\\f Ea \\o^ e + ||.9 Bc ||oa) !(</>,VOIk- 


Similarly, we have 



By Young’s inequality with constant /i e (1,2), we have 



5 e \v e - u e \ 2 dx > f i5 £ (|u £ | 2 - 2 |u e | |u £ | + |m £ | 2 ) dx 


>(l-M)lk £ iloA + ( 1 -^ 1 )l« e loA- 


Then, by the assumption specifically for (CDD), we have 6*3 > K, and so 



(4.6) 


+ (03 + K( 1 - v)W\\l Se + K( 1 - ^)\\u% Se 
>C(0i,K,ri ||(« e ,« e )|-t.- 


Hence, by the Lax-Milgram theorem, for every e > 0 there exists a unique pair of 
functions (u e ,v E ) e X E that is a weak solution to (CDD) and satisfies 



where the constant C is independent of e. Next, by the assumption that < 1 and 
Lemma 3.3, there exists a constant C, independent of e such that 


ll(«> E )lk < C(\\f Ea \\o^ + ||<Ao,0 <C(\\f Ea \\ LHn) + || 5 || i2(r) ), 


which is ( 2 . 22 ). 

By a similar argument involving the Lax-Milgram theorem, one can show the 
wcll-posedness of (SDD) in H 1 (fl,5 e ), (RDD) in V E , (DDDH) in W e and (NDDH) in 
iL 1 (H,^ e ). The independence of e for the constants in (2.23), (2.24), (2.25) and (2.26) 
follow again from Lemma 3.3 and the property that < 1. 

4.4. Weak convergence. Let ip € 7J 1 (H*) and ip 6 -n (T) be arbitrary. Let 
p Er e H 1 ( fl) denote the extension of ip to fl by the extension theorem [18, Theorem 
1, p. 254], and let 1 p Ec e 1L 1 (H) denote the extension of ip to fl as outlined in Corollary 


3.1. 


By < 1, Lemma 3.5, and Lemma 3.3, we see that ip Er e V e and ip Ec e H 1 (f2, S e ) 
for all e e (0,1]. Moreover, there exists a constant C > 0, independent of e, such that, 


b Er \\v e ^ C|Mliri(n*), |V ,£7c ||i,5 e < C\\ip\\ H i(ry 


Thus, we may test with tp Er and ip Ec in the weak formulation for (CDD). 

For e e (0,1], let (u e ,v £ ) e X e denote the unique weak solution to (CDD). Then, 
they satisfy 


(4.7) 


[&A Ea Vu e ■ Vp Er + &a Ea u E v Er + 5 E B Ec \7v e ■ \]ip Ec + S e b Ec v e ip Ec dx 

Jn 
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We analyse the bulk and surface terms separately. From (2.22), we have 

(4.8) flu e ||? iC . + ||u e ||g A < ||(</,« £ )ll* e < C(\\f Ea f L2(n) + ||<?||i 2(r) ) =: Cf, g , 

where C/, g is independent of e. Then, by Lemma 3.8 and Rcllich-Kondrachov com¬ 
pactness theorem, there exists u e H 1 (tt*) such that, along a subsequence 


(4.9) 

(4.10) 


u e |n* -- u in as e -* 0 , 

u e |n* -*■ u in L 2 (Q*) as e -*■ 0, 


(4.11) f 8 e u E (if,’ Ec - ip Er ) dx->• f 7 o(u)(^ - 7o(<^)) dU n 1 as e -* 0. 

Jn Jr 

By Holder’s inequality, f Ea ip Er e L 1 (H), and so by Lemma 3.4, we see that 

(4.12) f £ s f Ea ip Er dx = f £Jipdx + f £ e f Ea ip Er dx ^ f fipdx as e ^ 0. 

J n Jn* Jn\n * Jn* 

By Cauchy-Schwarz inequality, Lemma 3.4, (4.8), and the fact that (1 - £ e ) < \ < 
in Cl*. we have 


(4.13) 


fj 1 - «!“‘ IM dx 5 [L h"'! 2 dx ) : ‘ (L 11 ■ £l) ^ dx ) 1 

J. 1 

^ 1^1°,^ (^(!-^)b| 2 dx) <Cf, g [J^(l-£;e)\<p\ 2 dx) -»■ 0 as e -► 0. 

Thus, together with (4.10), we have 

I f £ e a £;o ii E <^ £> dx - f auipdx 
I Jn Jn* 

< f |a| \u e - €t| \ip\ dx + f |a| (1 - £ e ) |u E | |y>| dx + f la^l |u E | |</5 Er | dx 

Jn* Jn* Jn\n* 

< C a (iK - w|| i 2(n*)||<^||L 2 (n*) + 1 - £e) \u e \ \y\ dx + Cf,g\\<p Er ||z,2(n\n*,f e ) j 0, 

as e -*■ 0. This implies that 

(4.14) J~ t; e a Ea u e ip Er dx-» J~ auipdx. 

A similar calculation to (4.13) gives 

1 

J nt ( l-^ e ) |Vu £ ||V</?| dx < C ft g (j^l-^W dx j 0 as e - 0. 

Meanwhile, by Cauchy-Schwarz inequality, (4.8), and Lemma 3.4, 

I f &\7u s ■ (A Ea ) T \7ip Er dx < ||VM £ |U2 ( n.n*,c £ )C , (A Ba )||V^ Er |U 2 (nvn%c e ) 

I JS 2\f2* 

- Cf,gC(A Ea )\\ V^ Sr ||L 2 (nsn*,{ e ) -* 0 as e ->■ 0. 
Then, together with the weak convergence of Vu E |n* to Vu in L 2 {Cl*), we obtain 

I f £ e A Ea Vu £ ■ SJip Er dx - f AVu-V^dxl 
| J£} Jfl* 

(4.15) < \f V(u £ - u) ■ A T Vpdx + ||A||z,~(n*) f (1 - £ e ) |Vw £ | |V</?| dx 

|Jn* I Jn* 


If ^u s -(A Ea ) T ^ip Er dx 
I Jn\n* 


0 as e -*■ 0 . 
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Next, by (3.36), we have 

\fj E {\ g ^ + ri-\ 9 E i-\r c \ 2 )^ 

^ J^Gs + V'l 2 - \g\ 2 ~ \ip\ 2 )dH n ~ 1 = gipdH 71 ' 1 ase-^0 . 

\\v £ \\U<K^y)\\ 2 Xe <Cf, g , 

and so the conditions of Lemma 3.9 are satisfied. Thus, there exists a function v € 
iL 1 (T) such that 

f S e b Ec v s ip Ec dx ^ f bvipdH 71 - 1 , 

(4.17) Jn Jr 

/ 6 e B Ec yv e • y*p Ec dx / By tv • VrV’d'H™ -1 , 

Jn Jr 

as e -*• 0. Moreover, using the fact that ip Ec - (p Er e H 1 { f2), we obtain from (3.86)i 
with b Ec = 5= 1, 


f 5 s g Ec ip Ec dx = 
Jn 


(4.16) 

By (2.22), we have 


(4.18) [ 8 s v e (ip Ec - <p Er ) dx -*■ [ v(ip - 7o(y>)) dH" 1 as £ -»■ 0. 

Jn Jr 

In summary, from (4.11), (4.12), (4.14), (4.15), (4.16), (4.17), and (4.18), we see 
that passing to the limit e -* 0 in (4.7) leads to 

/ AVu ■ yip + auipdx + / Byrv ■ yrip + bvib dH 71 - 1 

Jn * Jr 

+ [ K(v-'y 0 (u))(i/.’-"/ 0 ( l p))dH n ~ 1 ~ [ f‘pdx-[(3gi/jd'H n ~ 1 =0. 

Jr Jn* Jr 


In particular (u,v) is a weak solution of (CSI). But by the well-posedness of (CSI), 
we must have that u = u and v = v, and that the whole sequence converge as e -*■ 0 . 

We remark on how to prove Theorem 2.3 for the other problems: 

• For (SDD), it follows directly from (4.16) and (4.17). 

• For (RDD), it follows from (4.12), (4.14), (4.15), (4.11) with ip Ec - ip Er re¬ 
placed with (p Er , and using Corollary 3.3 gives 

[ 5 e g Ec ip Er dx -»■ f gj 0 (ip) dT~L n ~ 1 as £ ->■ 0. 

Jn Jr 

• For (NDDH), it follows from (4.14), (4.15) and (4.12) with f Ea replaced with 
f Ea + V • (. A Er yh Er ) - a Ea h Er . 

• For (DDDH), we consider an arbitrary <f> e Hq{ il*) and extend it to <j> Er € 
iJp 0 (fi) by the reflection method. Then, by Lemma 3.7, <p Er e W e is a valid 
test function. From (2.25) and Corollary 3.2, we obtain, along a subsequence, 

w e D |n* wd e H 1 ^*) as e -» 0 , 


for some wd e Hq(Q*). Moreover, by Cauchy-Schwarz inequality, (2.25), and 
(3.50) applied to <f > Er , we see that 


I / 1 

|Jn e 


5 e w e D (/) E ' dx 


< C(\\f Ea \\ L 2 (n), \\g Er IIi/ 1 (O))II <p Er IIo,i, 5 e 0 as e 0. 


Then, together with (4.12), (4.14), (4.15), and the above, we see that u>d 
satisfies the weak formulation for (DSIH). 
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4.5. Strong convergence. We can choose ip - u e H 1 ^ 1 ) and ip = v e H 1 ( T), 
where ( u,v ) e x 77 1 (r) denote the unique weak solution to (CSI). Then, 

the extension u Er e by the extension theorem [18, Theorem 1, p. 254], and 

ip Ec e as in Corollary 3.1 are admissible test functions in the weak formulation 

of (CDD). Due to the coercivity of the bilinear form (4.4) (see (4.6)) and that | 
in Q*, we obtain 


lcDD«u e - u Er y - v Ec )) - a C DD((u Er , v Ec ), ( u e - u E \ v e - v Ec )) 
= ac DD {(u E - u E \v e - v Ec ),(u e - u Er ,v e - v Ec )) 


> C{6 U K) (| u e - u E r\\l ie + || U e - u E % Se +|K - 


> C(0,, JT) (K In. 


+ u - V 


Ec II2 


,a) • 


We claim that the left hand side converges to zero as £ -* 0. Indeed, by (4.12) with 
< p Er = u Er , (4.14) with a Ea = a = 1, and replaced with f Ea , we have 


(4.20) 


f £ >e f Ea (u e -u Er ) dx-*• f f(u - u) dx = 0 as e -* 0 . 
Jn Jn* 


Similarly, by (4.16) with ip Ec = v Ec and (3.85), we have 


(4.21) f 5 e g Ec (v e - v Ec ) dx f g(v - v) dJT " 1 = 0 as e - 0. 

Jn Jr 

Thus, 

lcDD((u e -u Er ,v e -v Er )) = f (ef Ea (u e - u Er ) + S e /3g Ec (v E -v Ec )dx -+ 0 as e 0. 

Jn 

Meanwhile, using (4.12) with f Ea tp Er replaced with a Ea \u Er \ 2 and with A Ea SIu Er ■ 
Vu Er , (4.14) with ip Er = u Er and (4.15) with cp Er = u Er , we have 

(4.22) f tp e {A Ea \Iu Er ■ V(u e - u Er ) + a Ea u Er (u e - u Er )) dx -*• 0 as e -»■ 0. 

Jn 

Next, by (3.35) with f Ec = v Ec , (3.86) with ip = tp Ec = u Bc , we have 


(4.23) 


[ S e (B Ec Vv Ec ■ V(v e - v Ec ) + b Ec v Ec {v e - v Ec )) dx - 0 as e - 0. 

Jn 


Finally, by (3.40) applied to / = u Er (2v Ec -u Er ) e H2 rl - 1 (Q), (3.54) with <p = v Ec -u Er e 
(3.35)i with = 1 and / Bc = (3.85) with and (3.86)i with 

b Ec = 1 and tp = we obtain 

f 8 Jv Ec - u Er ) (u e - u £ - + w Br ) dx 

Jn 

X c ( Er / o Ec Er\ e/ Ec Er\ I , e Ec e Er\ 

0 e (u (2v -u )-u (v ~U ) - |v \ +V V -V U J 
-»■ J t Jo(u)(2v - 7q(m)) -7o(m)(u-7o(m)) - |u] 2 + M‘ - «7o(w) dH" -1 = 0 as e -» 0. 


Thus, 

acDD (( u Er ,v Ec ) : {u e - u Er ,v e -v Ec )) 

= f t e (A Ea Vu Er ■ v(u e - u Er ) + a Ba w Br ( M e - u Er )) dx 

Jn 

+ J' 6 e (B Ec Vv Ec ■ V(v e - v Ec ) + b Ec v Ec (v e - v Ec ))dx 

+ K [ S e (v Ec - u Er ){v e -u e - v Ec + u Er ) dx - 0 as e - 0. 
Jn 
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From (4.19), this implies that 

(4.24) ||w e |n* -u\\ 2 H i(n*) + IK - u Er \\o,s e + K - ^lli ,s e ^ 0 as e ^ 0, 

which is the first assertion of Theorem 2.4 for (CDD). 

Furthermore, by the triangle inequality, (4.24), and (3.36) with f Ec = v Ec , we 
obtain 


| IK 111,<5* - IMI-fU(r)| ^ |IK|i,«. - Ill'll i,5 e | + |!l u ' Ec |i,«* “ IMK(r)| 

< II- v Ec \\i t s e + 11 v Ec 1 i ; 5 e - |n||jji(r)| -*■ 0 as e -»■ 0 . 

Hence, we obtain the normal convergence 

IKIIm, -► IMK(r) as e 0 . 

2 

Similarly, by (3.40) with / = \u Er \ , we have ||w Br ||o,< 5 c || 7 o(^)||L 2 (r) as e -> 0. 

Together with the triangle inequality and (4.24), we have 

Kllo,^ -* ll7oKK(r) as £ -* 0. 

We remark on how to prove Theorem 2.4 for the other problems: Recall that 
v E s e 5 e ), u R e V E , w e N e H 1 (H,^ E ) and w e D e W E denote the unique weak 

solutions to (SDD), (RDD), (NDDH) and (DDDH), respectively, and vs e i7 1 (r), 
ur,wn e F7 1 (H*) and wd e Hq(CI*) denote the corresponding unique weak solutions 
to (SSI), (RSI), (NSIH) and (DSIH). We consider the extensions v Ec as constructed in 
Corollary 3.1, u R r ,w Er by the method of reflection and w R ° € 0 (f2) obtained from 

extending wd by zero outside Cl*. Then, by Lemma 3.3, Lemma 3.5, Lemma 3.7, and 
that £ e < 1 in Cl, we have that v Ec e i7 1 (H,5 E ), u R r ,w Er e V E and w E0 e W E are valid 
test functions for the weak formulation of (SDD), (RDD), (NDDH) and (DDDH). 

• For (SDD), it follows from the above analysis by setting v e = v e s , v Ec = v Ec , 
P = 1 ,K = 0, A Ea = 0 nxr \ a Ea = 0 and f Ea = 0 . 

• For (RDD), it follows from setting u e = u e R , u Er = itf r , B Ec = 0 nxn , b Ec = 0, 
v E = v = 0, K = /3 and using (4.20) and (3.96) to see that 

Irdd(u e r - u R r ) := £,sf Ea (u e R - u R r ) + (35 E g Ec (u R - u R r ) dx 0 as £ -* 0. 

• For (NDDH), it follows from setting u e = w%, u Er = w% r , B Ec = 0 nxn , b Ec = 0, 
K = 0, g Ec = 0 and replacing f Ea with f Ea + V • ( A Er S7h Er ) - a Ea h Er . 

• For (DDDH), it follows from setting u e = w e D , u Er = w%°, B Ec = 0" xn , b Ec = 0, 
K = i, g Ec = 0, and we use (4.20) and (4.22) to deduce 

Idddh{w e d - wf, r ) 

:= f n Uf Ea (w e D - <1 ~ A Ea S/g Er • V(w e d - w Er )) dx 

- [ l; e a Ea g Er (w E D - w Er ) dx -* 0 as £ -* 0 , 

J n 

and use Cauchy-Schwarz inequality, (2.25), and (3.50) to deduce 
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5. Discussion. 

5.1. Alternate approximations for the Dirichlet problem. In the deriva¬ 
tion of (DDDH), we substituted 3 - However, an equally valid diffuse domain 
approximation to (DSIH) is (DDDH) m : 

-V • (ti e A E V» £ ) + 6 a E w e + -^~S e w e = £e/ E + V • (Z e A E Vg E ) - &a E g E in 12, 

£ m 

A E SJw e ■ v = 0 on <912, 

where m > 0. We see that (DDDH) is the case m = 1. Let us introduce the weighted 
Sobolev space 


W™ := {/ : 12 -»■ R measurable s.t. / |n E e W l ' 2 { 12 e ,£ E )} with 

(/, 5 >w™ : = f &(/<? + V/- Vg)+ ~^6 e fgdx. 

e Jn e m 

One can use the Lax-Milgram theorem to show that (DDDH) m is well-posed in W™ 
with unique weak solution w e D m . 

For g e C 1 (fl) such that g |r= 0, we see that for any m e (0,2), by (2.10) with 
q = 1 + to, and (3.42) with q = oo, 

f — S e \g\ 2 dx ->• 0 as e ->• 0 . 

Jn e m 11 

However, due to (3.42) with q = 2, for any / e H p 0 (D), we can only show that there 
exists a constant C > 0, independent of / and e € (0, 1 ], such that 

f Q I/I 2 dx ^ c \\f\\ 2 H^n) for m 6 (°- !]■ 

Thus, for any m e (0,1], using the techniques presented here, we can show 
\\ w D,m In* ~ w d ||jyi(fi*) ^ 0 as £ ^ 0, 

where W£> is the unique weak solution to (DSIH). We point out that it is not known 
if the formal asymptotic methods used in [1, 22, 31] can be applied to non-integer 
to. However, this result may be beneficial for numerical computation, where high 
values of to and small values of £ usually result in stiff computations, and thus are 
computationally expensive. 

5.2. Issue of uniform estimates for (NDD). Let us define 

dNDD(u,4>) : = InDd(<P ) : = l e B(f Ea A) + i%(.g Ec A >)■ 

Suppose we seek weak solutions in H 1 ( 12,£ E ), then o,ndd (•, •) is bounded and coercive 
uniformly in e. But, by (3.53), the boundedness of Indd over 7L 1 (f2,£ e ) scales with 
j. Thus, we cannot obtain uniform estimate similar to (2.26). If we seek solutions in 
V e , then Indd is bounded uniformly in £. But cindd is not coercive in V e uniformly 
in £. Hence, to obtain uniform estimates in £, we have to remove the term involving 
5 e , which amounts to transforming (NSI) to (NSIH) through Theorem B.2. This, 
then motivates Assumption 2.4 and the corresponding diffuse domain approximation 
(NDDH) now has a weak solution in iL 1 (D,^ £ ) with uniform estimates in £. 
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We point out that if a trace-type inequality holds: 

f 5 e \uf dx < C f £ e (|u | 2 + |Vu| 2 )dx, 

Jn Jo, 

for some C > 0 independent of e, then the Hilbert spaces V e and are 

equivalent uniformly in e. Moreover, (NDD) will be well-posed over H 1 ( 12, £ e ) = V e , 
and Assumption 2.4 is not required. However, stronger assumptions than Assumptions 
2.7 and 2.8 are required, see for instance [11, Theorem 4.2] or Section 5.5. 

5.3. Regularity of T for (CSI) and (SSI). In the proof of Lemma 3.9, we 
used that the components g l J , 1 < i,j < n, of the inverse of the metric tensor Q~ 3 are 
C 1 functions. Then, from the definition, we require that T is of class C 3 . Since we 
apply Lemma 3.9 only to (CDD) and (SDD), for the other problems, a C 2 regularity 
for T is sufficient. 

We point out that using (3.26) and (3.14), we can consider the following splitting 
in (3.88): 


n — 1 

(Vr-F £ )r= 'Zsi r da i V e d ar G e 

(5-1) J= \ 

= Yj g° Q r d Sj V e d Sr a + ezgl r d Sj V e d Sr a + £jrd Sj V e ds r Ge. 

3 =1 

If we have some control over the components {d Sj V e )i<j< n -i, then it is sufficient 
to apply integration by parts only on the first term. This in turn implies that we can 
potentially drop the required regularity of T from C 3 to C 2 . However, the hypothesis 
that ||v e ||i ) 5 e is bounded uniformly in £ seems to be not sufficient to give any control 
over the components (d Sj V e )i<j< n -i- 

5.4. Comparison with the results in [16]. In a time-independent setting, 
by choosing q - 1, p e (s,z ) = ~p{z) = <5(z), the function spaces B and X employed in 
[16] are equivalent to L 2 (Tub £ (T), 5 e ) and R 1 (Tub e (r), S e ), respectively. Moreover, 
by comparing with the notation and results in Section 3.2, the results of [16] in our 
notation are 


c e (s, z) c(s) e B <=> |c £ - c|| i 2 ( Xub e( r ) -»■ 0, 


and 

c E (s, z) ^ c(s) e X / S e (c e - c)tl> dx -»■ 0, / S £ V(c e - c) • Vifidx -»■ 0, 
JTub E (r) JTub e (r) 

for if) e 7L 1 (Tub £ (r),<5 e ). Thus, in the limit £ -> 0, c e converges weakly to a function 
c defined only on the surface T. 

We point out that, in the proof of Lemma 3.9, it is crucial that (3.90) holds, i.e., 
the test function ip is extended constantly in the normal direction. Otherwise, for an 
arbitrary test function A e H ^If,^) with representation A e in the (p, z) coordinate 
system, a similar calculation to (3.87) yields that 

X(-ri e , k ,ri eik )(z)-d z A E is bounded in L 2 (T x R,J) for all £ s ( 0 , 1 ]. 
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When computing (3.91), we have an additional term of the form 

(5.2) f f 6(z)^d z V £ d z A £ \l+£zC H (p) + C R (p,£z)\ d'H” -1 dz, 

Jve.k J r £ z 

and by (3.87) and Cauchy-Schwarz inequality, we obtain only the uniform bounded¬ 
ness of (5.2) in £, and we are unable to show that (5.2) converges to zero as £ -* 0. 

Thus, we are unable to enrich the space of test functions of the weak convergence 
in Theorem 2.3 and Lemma 3.9 from Cp (2.27) to i7 1 (S2,<5 e ). However, in [16, Proof 
of Theorem 4.1], to show that the limit c satisfies the correct weak formulation of an 
advection diffusion surface PDE, the authors considered a test function x e X with 
d z x = 0, which is similar to what we do in Lemma 3.9. 

5.5. Comparison with the results in [11]. The framework developed in [11] 
bears the most resemblance to the framework presented here. In particular, weighted 
Sobolev spaces are employed for the well-posedness of the diffuse domain approxima¬ 
tion, and the authors also studied tubular neighbourhoods and derived results similar 
to those in Section 3.1. In our notation, the object of study in [11] is the Robin 
boundary value problem 

-V • (AX/u) + cu = / in fU, 

5.3 , ' 

AX7u • v + bu = g on T. 


With the choice S e = |v£ e |, the diffuse domain approximation in weak formulation is 
given as 


(5.4) 



u e ■ Vv + cu e v) + ]V£ e | bu e v dx 



gv dx 


for all v e 7L 1 (fl e ,^ e ), where f l e is defined in (2.13). Moreover, the authors can 
show a trace-type inequality: There exists a C > 0 such that for all v e W 1 ’ p (Q £ ,f; E ), 
1 < p < oo, and £ sufficiently small, 

(5.5) f M P |V&|d x<C f £s(M P + |Vvndx. 

Jtt e J 

Additionally, it is assumed that satisfies the following behaviour: There exists 
Ci, £2 >0 and a > 0 such that for all x e f2 e , 


(5.6) 


/dist(a:,9H e ) ’ 




Then, the authors can show continuous and compact embeddings from 1 / P 1 ’ p (f2 e ,£ £ ) 
into L^n^Ce), as well as a Poincare type inequality: There exists a Cp > 0 such that 
for all v e W l ’ p (Xl £l £e), 1 <p < 00 and e sufficiently small, 


(5.7) 


p 


< Cp ( |Vti£ P(fWe) + J ^ |tf |V6I dx ). 


Using these results, together with the fact that f l £ = f2*uT e and |P e | < Ce, the authors 
deduced that 


/ dx - / h dx 

Jn, Jn* 


C'l^llip(r e ,5 e )£ rl p , 

C\h\ w i,p(r e ^ e )£ t 
-> 0 as £ -»• 0, 


for h e L p (fU,£ e ), 1 <p < 00 , 
for h g 1 <p < 00 , 

for h e L 1 (U). 
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We note that the last case is similar to Lemma 3.4, and both proofs use the dominated 
convergence theorem and the pointwise limit of £ E . 

By the application of (5.5), (5.7) and the Lax-Milgram theorem, the diffuse do¬ 
main approximation is well-posed with a unique weak solution u E € iJ 1 (f2 E ,£ E ). The 
chief results of [11] is the following error estimate: 

(5.8) ||u-u £ ||wi.2(fi e ,j e ) < Ce^-p, 

where u e W 1,P (S2*), 2 < p < 2* is the weak solution to (5.3), u e e is the 

unique weak solution to the diffuse domain approximation, C is independent of e, a 
is the exponent in (5.6), and 


P 


* 

a 


p[n + a) 
n + a- p 


if p < n + a, 


Pa~°° if p>n + a. 


We point out the differences between our work and the work of [11]: 

• [11] utilises a double-obstacle regularisation and the fact that |T E | < Ce for 
some constant C to deduce (5.8). In our work, we cover both the double-well 
and double obstacle regularisations. 

• We do not have a rate of convergence for the error \\u e R - Ujf r ||i,£ e , but [11] 
requires W 1,p regularity of the original problem to deduce (5.8). This follows 
in the same spirit as (3.41), where for more regular functions, we are able to 
deduce a rate of convergence. 

• For the choice S E = |V£ e |, the trace-type inequality (5.5) allows [11] to handle 
non-homogeneous Neumann boundary problems. Our setting is more general, 
where S E and £ E can be unrelated, as long as (2.11) is satisfied. But as 
discussed in Section 5.2, we require stronger assumptions on the data to 
compensate for the generality. 

• In [11], the non-homogeneous Dirichlet problem is handled with a penalization 
method in the same spirit as Lemma 2.1. The corresponding diffuse domain 
approximation is (5.4) but with j4r|V£ e | replacing |V£ e | with m > 0. The 
error estimate for the Dirichlet problem becomes 

IIIn* ~ u d\\h 1 (q*) ^ C{e m + e 2 p m ), 

where u e W 1,P {VL*), 2 < p < 2*, is the unique solution to (DSI), and 
u e D e H 1 (f2 e ,£ e ) is the unique solution to the diffuse domain approxima¬ 
tion. An optimal choice is m = i - leading to a rate of convergence of 

order C?(e 3_ ^). In contrast, we homogenise the problem and use Lemma 2.1 
to derive a family of approximation (DDDH) m , m > 0. Although we do not 
have any rate of convergence, we showed strong convergence for any m e (0,1] 
(see Section 5.1). 

• Lastly, in this work, we are also able to consider partial differential equations 
purely on the boundary T or coupled with the bulk equation, whereas [11] 
only focus on a bulk elliptic equation with Robin, Neumann and Dirichlet 
boundary conditions. 

Appendix A. Derivation of the distributional forms. A key result we use 
to derive an equivalent distributional form is that of Alt [5, §2.7 and Theorem 2.8] 
(see also [22, Theorem 6.1]). We shall give a time independent version: 
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Theorem A.l. Given an open set n c R n containing an open subset Cl* such 
that T = dD* c Cl. For itT, let v(x) e (X^X ) 1 be the external unit normal of I\ 
Denote by xn* and dr the following distributions: 

f f d X n*= f f dx, f fdS r = f /dH”' 1 . 

JQ Jq* Jo, Jr 

Then, for smooth functions q h j, f b : SI* -» M, gj, / s : T -> R ; 1 < j < n, the distributional 
law, 


V- (xsW' + fW 5 ) =Xn*/ b + <5r/ s inV'{Cl), 


is equivalent to the following: 

V-q b = f b inD*, 

Q. S {p)^T p Y, S] T -q s = f s + q b -v onT, 

where Vr( - ) denotes the tangential gradient on T. 

To derive the equivalent distributional form for (CSI), we set 

q b = -AS7u, q s = 0, f b = / - au, f s =K(v-u). 

Then, by Theorem A.l, 

-V • (xn *A\/u) + xn *au = Xn*/ + S r K(v - u) in V'{Cl) 


is the equivalent distributional form for 

-V • (AVu) + au = / in fX, 

AVu ■ v = K(y - u) on I\ 

Now, setting 

q b - 0 , q s = -B\7v e , f b = 0, f s =/3g-bv-K(v-u), 


where v e is the extension of v constantly in the normal direction as outlined in Section 
3.3 (see also (3.32)). Then, 

-V • (dr£>Vu e ) + Srbv = Srfdg - SrK(v - u ) in V' (Cl) 


is the equivalent distribution form for 

-Vr • (BVrv) + bv + K(v - u) = /3g on F. 

Hence, the equivalent distributional form for (CSI) is 

-V • (xn *A\7u) + xn *au = xn */ + S r K(v - u ) in V'(n), 
-V • (dr BSJv e ) + 6rbv = 6r(3g - SrK(v - u ) in V'{n). 


Appendix B. Properties of the trace operator. We use the following result 
from [25, Theorem 1.5.1.2] regarding the continuous right inverse of the trace operator 
for Sobolev functions (see also [43, Theorem, p. 211] and [35, Theorem 3.37]): 
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Theorem B.l. Let k > 0 be an integer and 0* is a bounded open subset of 1" 
with C k ^ boundary T. Assume that s - i > 0 is not an integer with s < k + 1. Then 

the trace operator 70 : W s ’ p (fT) -*■ W S ~ 2 ’ P (T) has a continuous right inverse. 

The derivative AVu ■ v on T is defined by 

n 

AVu ■ v := ^ Vj^ 0 (AX/u)j on T, 

3 =1 

where v is the outward unit normal to I\ This object is well-defined for u e H 2 (Tl*) 
and A e (C°(fl*)) nxn . The analogous result to Theorem B.l for the derivative AVu-u, 
is a consequence of the solvability and regularity of solutions to non-homogeneous 
Neumann problems (see [35, Theorem 4.18 (ii)] and [40, Theorem 4]). We also refer 
the reader to [25, Theorem 1.6.1.3] for a similar result with stronger assumptions on 
A. 

Theorem B.2. Let k > 0 be an integer and f 1* c R ra is a bounded open set, 
with boundary T e C fc+1,:L and outward unit normal v. Assume that the matrix A(x) 
is uniformly elliptic and each coefficient belongs to C k,1 (tt*). Then, if f e 
g e H k+ z(T), the non-homogeneous Neumann problem (with X > 0) 

-V • (AVu) + A u = f in f2*, 

AlVu • v - g onT, 

admits a unique solution u e H k+2 (Tt.*) with 

^ CGMI-fUCfi*) + ll/llH fc (n*) + |p| 

We will use the following compactness property of the trace operator 7 o(-) from 
[4, Theorem A 6.13, p. 257] (see also [14, Theorem 3.8.5, p. 167] or [38, Theorem 6.2, 
p. 103]): 

Theorem B.3. Let Q c R™ be an open bounded domain with Lipschitz boundary, 
and 1 < p < 00 . For k € N, let Uk,u € W 1,p {Tl). Then, as k -* 00 , 

Uk -- u in TT 1 ’ p (fl) => 7o(rifc) -*■ 7o( u ) in L p (dfl). 
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